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Algoritmes  voor  het  fuseren  van  twee  sets 
(sonar)  data 


Bij  het  combineren  van  meerdere  sensoren  is  het  belangrijk 
de  informatie  van  de  verschillende  sensoren  20  goed 
mogelijk  samen  te  voegen.  Dit  kan  grote  impact  hebben  op 
de  prestatie  van  het  gecombineerde  systeem.  In  dit  rapport 
worden  verschillende  manieren  onderzocht  om  de  informatie 
van  de  verschillende  sensoren  (in  ons  geval  sonars)  te 
combineren. 


mR{b(ue)P  us  {red),  (black) 
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P  rob  leemstel  ling 

Een  groat  dee  I  van  het  programma 
* Onde  r/ee  waterbee  Ido  pbou  w  in  ondiep 
water'  (V506)  is  gericht  op  ondeczoek  op 
hei  gcbicd  van  multistaiische  operatic  s. 

De  Koninklijke  Marine  heeft  hier  inieressc 
in  omdai  dcrgelijke  operaties  komende 
jaren  /owe  I  naiionaal  aLs  internal  ionaal  een 
Steeds  belangrijkere  ml  zullen  gaan  spelen 
t  ijde  ns  onde  r/eebi  x  nhest  rijd  i  ngsi  ipe  rai  ie  s, 
Bij  multi  statisehe  operaties  worden  sonars 


van  meerdere  platformen  tegelijk  ingezet. 
wat  grote  tactisehe  voordelen  kan  hebben  en 
een  verbeterde  prestatie  be  loo  ft.  Voor  een 
op ti male  prestatie  in  een  multistatische 
operatic  diem  de  verk regen  informal  ie  van 
alle  sonars  naar  een  eemraal  punt  gestuurd 
te  worden.  waar  het  gefuseerd  kan  worden. 
In  dit  rapport  wordi  hei  pmhleem  van  de 
fusie  van  (sonar)data  be  s  pm  ken.  in  de 
aanwezigheid  van  fouten  op  de  po&tties  van 
getde nt i ficee rde  (sonar )co nt act e n . 


De  comacien  zijn  over  het  algemeen  van 
verschillende  sonarsystemen  afkomstig  en 
worden  meestal  niet  op  p redes  de  goede 
fxssitie  waarge  nomen.  Hoe  kunnen  de/e 
opiimaal  worden  gecombineerd  om  een  m 
gmoi  mogelijke  prestatie  of 
p  re  static  ve  rbetering  ten  opzicht  van  de 
i  nd  i  v  id ue  le  so  nars  te  he  h  he  n  ? 


Kesultatcn  en  concludes 

Uit  de  theorie  volgl  dat  de  posit  te  van  een 
contact  na  fusie  significant  anders  kan  /ijn 
dan  de  posit  ie  van  de  oorspronkelijke 
contacten  waarnit  hei  is  vtxmgekomen. 
Deze  positie  na  fusie  ligt,  tegen  de 
verwaehting  in.  meestal  niet  op  de 
verbindingslijn  lussen  de  twee 
0<> rsp  mnke  1  ijke  objec  te  n 
De  deteclie prestatie.  na  satnenvoegen  van 
de  doeleontaeten  en  loosalarmen  van  twee 
sonars,  verbetert  aanzienlijk.  met  I  tot  5  dIV 
Daadwerkelijk  fuseren  van  de  twee  sets 
contacten  (dat  wil  zeggen.  in  de 
samen  go  voegde  set  de  contacten  die  ‘die  hi 
genoeg  bij  elkaar  Liggen’  combineren  in 
nieuw  contact)  voegt  hier  geen  verbetering 
aan  toe  in  terrnen  van  detectiepresiatie 
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Enkele  van  de  onderzochte  fusie-algoritmes 
laten  echter  wel  een  verheierde 
positiebepaling  van  de  gefuseerde  contacien 
zien  ten  opzichie  van  die  van  de  originde 
contacted  Dit  kan  belangnjk  ztjn  als  de/.e 
contacien  bijvoorbeeld  aan  een 
trackingaJgoritme  worden  gevoed. 
aangezien  do  prestatic  van  een  dergelijk 
algoriime  sterk  afhangt  van  nauwkeurige 
posit ie-informatie  van  dc  invocr. 
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Algontmes  voor  het  fuseren  van  twee  sets  (sonar)  data 


Toepasbaarheiri 

[Je  bcschrcvcn  fusie-algoritmes  kunnen  in 
dc  hestaande  (IR)LFAS- 
signaalverwerki ngsketen  ingebouwd 
worden.  Als  cr  dan  muJtistattsche  operatics 
worden  uitgevoerd  met  het  (IR)LFAS- 
systeem  kunnen  de  algontmes  een  gewenste 
tussenstap  vormen  aivorens  a  lie 
sonafcomacten  aan  een  trackings  Igorit  me  te 
voeren. 


PROJECT 

Projectbegeleider 
ing.  E.  van  de  spek, 
DMQ/DWS&B/RZS&B 

Projectl  eider 

dr.  P.  A.M,  de  Theije,  TNO  Defensie  en 
Veiligheid 

Projecttftel 

LFAS  netcentric 


Contact  en  rapportinformatie 


Oude  Waalsdorperweg  63 
Postbus  96864 
2509  JG  Den  Haag 

T  +31  70  374  00  00 
F  +31  70  328  09  61 

info-DenV@tno.nl 

TNO-rapportnummer 
TNQ-DV  2006  A518 

Opdrachtnummer 


Datum 

december  2006 

Auteur(s) 

dr.  P.A.M,  de  Theije 
if.  C  A  M.  van  Moll 

Rubricering  rapport 
Ongerubriceerd 


Projectnummer 

015.34675 

Projectplanning 
Start  01-01-2005 
Gereed  31-12-2008 


Projectteam 

dr.  P.A.M.  de  Theije, 
ir.  C.A.M.  van  Moll 


ONGERUBRICEERD 


ONGERUBRICEERD  |  TNO  report  |  TNO-DV  2006  A518 


4/51 


Abbreviations 


JRP 

Joint  Research  Project 

ML 

Maximum  Likelihood  (estimate) 

NURC 

NATO  Undersea  Research  Centre 

pdf 

probability  density  function 

ROC 

Receiver  Operating  Curves 

SNR 

SignaMo-Noise  Ratio 

TNO 

Netherlands  Organisation  for  Applied  Scientific  Research 

TNO-DenV 

TNO  Defence.  Security  and  Safety 
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List  of  symbols 


Theory,  general 


N(p.<r) 

Normal  distribution  with  mean  fi  and  variance  cr 

R 

Receiver  system  1 

S 

Receiver  system  2 

vAj} 

Probability  of  event  A  /  or  A2 

p{Al  aA2\b} 

Probability  of  event  Ai  and  A  2  under  the  condition  of  event 

P{a  £  A} 

Probability  of  event  that  a  is  not  an  element  of  the  set  A 

£{.v} 

Expected  value  of  the  statistic  x 

Theory,  one-dimensional  case 

T 

True  target  position 

T< 

Stochastic  variable  of  target  position 

% 

ML  estimate  of  target  position 

WR 

Stochastic  variable  of  position  of  detection  by  receiver  R 

(Or 

Observation  of  Wr 

(Jr 

Standard  deviation  of  Wr 

w* 

Stochastic  variable  of  position  of  detection  by  receiver  S 

m 

Observation  of  Wx 

<?s 

Standard  deviation  of  Ws 

Theory,  two-dimensional  case 

T  =  (Xt ,  }t) 

True  target  position 

Tc 

Stochastic  variable  of  target  position 

%  =  (  TX  ,  TV) 

ML  estimate  of  target  position 

WK 

Stochastic  variable  of  position  of  detection  by  receiver  R 

<J>R  =  (-T/f .  Vft) 

Observation  of  Wr 

Cr 

Covariance  matrix  of  Wr 

ws 

Stochastic  variable  of  position  of  detection  by  receiver  S 

=  (.Vs .  Vs) 

Observation  of  W* 

cs 

Covariance  matrix  of  WH 
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List  of  symbols  (continued) 

Simulations 


c 

Sound  speed  in  water  [m/s] 

(•t.t.i  •  V.t.i) 

{A.r,, .  Ay,,,, .  an. 
A  0, 

A  (p, 

A  ti 

(x  ,  y  imposition  sonar  i  (i  =  1,2)  [m] 

Error  in  (x  ,  y)- position  sonar  i  [m] 

Error  in  receiver  heading  sonar  i  [deg] 

Error  in  bearing  estimate  with  sonar  i  [deg] 
Time-synchronisation  accuracy  sonar  i  [s] 

SNRj:l 

Ur.; .  y,(i) 

Target  signal- to- noise  ratio  of  target  j  relative  to  sonar  i  [dB] 
(x  ,  y  )- position  of  target  as  observed  by  sonar  i  [m] 

At 

AAt 

Amplitude  threshold  [dB] 

Increase  in  amplitude  threshold  [dB]  due  to  data  fusion 

NtUpiAr) 

N/a 

pdJAr) 

PfaAAr) 

Number  of  detections  by  sonar  i  for  ping  p 

Number  of  false  alarms  by  sonar  /'  for  ping  /; 

Probability  of  detection  by  sonar  i 

Probability  of  false  alarm  by  sonar  i 

Q 

Set  of  threshold  crossings  from  sonar  i 

ONGERUBRICEERD 


ONGERUBRICEERD  |  TNO  report  |  TNO-DV  2006  A518 


7/51 


Summary 


In  this  report  we  have  focussed  on  the  topic  of  multi  static  sensor  and  data  fusion. 

First,  we  have  established  a  theoretical  foundation  for  the  probability  that  two  contacts, 
as  observed  by  two  separate  systems,  are  due  to  the  same  target. 

This  derivation  includes  the  position  errors  in  the  observations  that  are  expected  for 
both  systems.  Apart  from  the  association  probability,  the  most-likely  target  position  is 
obtained.  The  theory  is  derived  for  the  one-  and  two-dimensional  cases. 


To  test  the  above  theory,  and  to  assess  die  benefit  of  data  fusion,  realistic  simulations 
have  been  carried  out.  These  simulations  are  based  mainly  on  common  Rayleigh 
statistics,  both  for  the  background  and  the  target  amplitudes.  They  take  into  account 
position  errors  in  the  observations  (to  mimic  the  fact  that  any  sensor  will  not  yield 
infinitely  accurate  contact  positions),  finite  sampling  (in  time  and  bearing)  of  the 
sonars,  and  can  include  targets  of  different  signal-to-noise  ratios  and  at  different 
positions. 


The  following  definitions  are  used  throughout  the  report: 

*  ‘Combining’;  All  contacts  identified  by  the  two  separate  systems  are  stored  in  a 

single  set  of  contacts, 

*  ‘Fusing*:  As  in  'combining’,  but  for  each  pair  of  contacts  from  both  systems 

it  is  checked  if  they  originate  from  the  same  target.  If  so,  this  pair 
of  contacts  is  merged  into  a  single  contact. 

*  ‘OR’  fusing:  All  contacts  that  have  formed  after  ‘fusing’  are  taken  into  account. 

*  '  AND’  fusing:  Only  those  contacts  that  have  formed  after  ‘fusing’  and  that 

originate  from  two  contacts,  one  from  each  system,  are  taken  into 
account. 


The  main  conclusions  from  the  simulations  are  as  follows: 

*  Combining  two  data  sets  of  independent  sonars,  the  total  detection  performance  of 
the  system  improves  significantly  in  comparison  to  a  single  sonar.  This  is  due  to  the 
enhanced  probability  of  detection,  despite  the  increase  (by  a  factor  of  two)  of  the 
number  of  false  alarms. 

*  Fusing  the  two  data  sets  with  the  'OR'  rule,  the  total  detection  performance 
(as  evaluated  from  the  ROC  curves)  of  the  system  is  equal  to  the  (improved) 
performance  of  the  combined  system.  Although  it  was  expected  that  fusing  two  data 
sets  would  give  a  better  performance  than  combining  them,  this  is  not  the  case. 

•  Fusing  the  turo  data  sets  with  the  'AND'  rule,  the  total  detection  performance  of  the 
system  decreases  relative  to  that  of  the  single  sonars.  The  reason  for  this  decrease  is 
the  inaccurate  positions  of  the  original  contacts,  wrhich  cause  only  few  pairs  of 
target-originated  contacts  to  actually  fuse.  Without  position  errors  included,  it  has 
been  shown  by  others  that  the  performance  after  ‘AND’  fusion  would  be  equal  to 
that  of  the  original  sonars, 

•  After  'OR'  fusion,  the  standard  deviations  of  the  observed  contact-position  errors 
decrease  with  threshold  Aj>  especially  at  higher  thresholds,  Ar>  16  dB. 

The  standard  deviations  have  roughly  reduced  to  one  third  of  their  original  value 
between  Aj  =  10  and  Aj  =  20  dB.  This  is  because  at  higher  threshold  relatively  more 
fused  contacts  are  present,  that  have  a  higher  position  accuracy. 
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•  After  'AND"  fusion*  the  standard  deviations  of  the  observed  contact-position  errors 
also  decreases  with  threshold  At,  die  decrease  being  more  profound  at  lower  values 
of  Aj  than  after  'OR'  fusion* 

The  following  conclusions  from  the  simulations  hold  for  all  fusion  algorithms: 

•  Applying  a  small  increase  in  the  threshold  A At  as  part  of  the  fusion  algorithms 
results  in  a  slightly  lower  ROC  curve  and  is  therefore  not  a  good  idea. 

•  The  performance  results  of  the  different  fusion  algorithms  do  not  change  noticeably 
with  the  number  of  targets  in  the  simulations*  at  least  in  the  range  1  -  30. 

•  For  high  mean  target  SNR  20  dB)?  some  fusion  algorithms  show  relatively  high 
false-alarm  probabilities  at  high  amplitude  thresholds  (At  >  16  dB)*  due  to  the 
fusion  of  target-originated  contacts  and  false  contacts. 

•  A  single  simulation,  where  the  target-amplitude  follows  a  one-dominant-plus- 
Rayleigh  distribution  rather  than  a  Rayleigh  distribution,  suggests  that  the  results 
quoted  above  do  not  depend  very  much  on  the  target-amplitude  distribution. 

The  above  results  indicate  that  combining  or  ‘OR* -fusing  two  datasets  gives  a 
considerable  increase  in  performances  in  comparison  to  a  single  sonar* 

However,  'AND'  fusion  is  not  a  good  option  in  terms  of  detection  performance* 
especially  not  in  the  case  of  large  contact-position  errors.  ‘OR'  fusion  does  not  add 
detection  performance  (compared  to  just  combining  the  two  data  sets)*  but  it  does  yield 
more  accurate  position  information  to  the  fused  data  set.  This  can  be  important  if  the 
resulting  data  sets  are  fed  to  a  tracking  algorithm*  especially  in  heavy  duster  areas* 
where  many  contacts  may  occur  in  a  small  volume. 
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1  In  trod  uci  ion 


LI  Background 

Multistatic  operations  can  have  great  advantages:  they  can  lead  to  larger  detection 
ranges  and  a  better  coverage  of  a  possible  submarine  operations  area.  Also,  such 
operations  can  be  less  dangerous  for  the  participating  units  as  only  one  source  is  needed 
and  all  other  units  can  operate  passively  and  thus  remain  covert.  Apart  from  tactical 
advantages,  also  acoustical  advantages  may  be  obtained  (e,g,  less  reverberation). 

In  addition,  significant  gain  is  expected  from  the  combination  of  mono-  and  bistatic 
detections.  Combining  the  information  of  multiple 

(either  monostatic  or  bistatic)  sonar  systems  should  in  theory  provide  enhanced 
detection  performance,  more  clues  towards  classification  [2],  and  increased  localization 
accuracy.  For  the  latter,  positioning  issues  play  an  important  role. 

In  [1  ]  [3]  [4]  it  was  shown  that  the  uncertainty  of  the  position  of  a  detection  by  a 
receiver  depends  in  general  on  die  errors  in: 

•  Sound  speed. 

•  Receiver  heading, 

•  Bearing  estimate. 

•  Accuracy  of  time  of  reception. 

•  Source-  and  receiver- position  errors. 

The  combined  effect  of  all  these  components  can,  in  some  cases,  generate  quite  large 
uncertainties  in  the  target-contact  position,  up  to  several  hundred  metres. 

The  relative  importance  of  these  parameters  depends  on  the  target  position  itself,  and  on 
the  receiver  being  monostatic  or  bistatic. 

For  an  overview  of  the  current  status  of  multistatic  research,  see  the  excellent  overview 
in  [5]. 

L2  Motivation  for  current  research 

The  current  research  is  carried  out  in  the  context  of  the  project  4  FAS  Netcentric  \  part 
of  the  programme  'Onderzeebootbestrijdingscapaciteii  in  ondiep  water*. 

The  main  themes  of  this  project  are: 

1  Data  fusion. 

2  Multi -sensor  target  cracking, 

3  Including  submarine  sensors  in  multistatic  operations. 

4  Netcentric  underwater  communication. 

5  Including  dipping  sonars  and  sonobuoys  in  multistatic  operations. 

This  report  focuses  on  the  first  theme:  Data  fusion. 

Trying  to  fuse  two  data  sets,  originating  from  different  receivers,  one  has  to  deal  with 
the  position  errors  as  target  contacts  are  in  general  not  perfectly  aligned. 

The  question  that  arises  then  is:  How  to  combine  in  an  optimal  way  the  contacts  from 
the  different  receivers ?  This  question  is  relevant,  e.g.,  for  multi  static  tracking 
algorithms:  if  multiple  sonar  systems  are  executing  a  joint  operation,  their  identified 
sonar  contacts  can  be  combined  and  fed  to  a  so-called  centralised  tracker. 

This  tracker  receives  all  sonar  contacts  of  all  sonars  and  from  these  tries  to  form 
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consistent  tracks  by  associating  contacts  at  different  time  stamps.  For  such  algorithms 
the  quality  of  the  input  is  very  important:  the  position  accuracy  of  the  input  contacts 
should  be  as  good  as  possible,  and  the  number  of  input  contacts  should  not  be  too  large. 
For  both  purposes  it  is  important  to  combine  the  contacts  of  the  different  sonars 
optimally. 

1.3  Outline  of  report 

This  report  is  organised  as  follows:  Chapter  2  describes  the  theory  of  how  to  associate 
two  contacts,  as  observed  by  two  sonar  receivers.  This  gives  an  association  probability 
for  these  contacts  to  be  due  to  the  same  target,  as  well  as  a  most-likely  position  of  the 
combined  contact.  In  Chapter  3  this  theory  is  illustrated  using  a  simple  example. 
Chapter  4  presents  the  simulations  that  are  used  to  assess  the  performance  of  the  data- 
association  algorithms,  as  well  as  an  evaluation  of  the  improvements  to  be  expected  due 
to  data  association. 
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2  Theoretical  analysis 


2-1  Statistical  concepts 

Consider  a  sonar  system  with  two  independent  observation  systems:  one  receiver  S,  and 
a  second  receiver  R.  For  example,  for  a  simple  bistatic  system,  S  can  be  collocated  with 
the  source  position,  while  R  is  assumed  to  be  at  a  certain  distance  from  the  source. 
However,  the  exact  operating  mode  of  S  and  R  does  not  change  the  algorithms  below-. 

A  signal  transmitted  by  the  source  reflects  on  the  target  and  may  cause  a  contact  by 
each  of  the  two  independent  observation  systems.  The  position  Ws  of  the  contact  by 
receiver  S  is  assumed  to  be  a  normally  distributed  stochastic  variable,  centred  at  the 
target  position  T  with  its  (co)variance  C$  determined  by  this  position.  The  position  \VR 
of  the  contact  by  R  has  its  owrn  normal  distribution,  also  centred  at  the  target  position 
hut  with  (co)variance  Cr.  That  is, 

Ws  -  Af(T,Cs)  and  WR  -  /V(T,CR) 

with  Mp  -  C)  representing  a  normal  distribution  with  mean  //and  (co)variance  C 
Since  both  contacts  are  assumed  independent,  the  probability  density  function  (pdf)  for 
a  combination  of  both  contacts  is  the  product  of  their  individual  probability  density 
functions.  However,  given  the  event  of  combined  detected  positions  (o)s  .  wr)  which 
come  from  the  same  target,  this  target  position  T  is  unknown  and  must  be  estimated 
from  this  given  combination.  This  can  be  realised  with  the  concept  of  the  maximum 
likelihood  (ML)  estimator,  by  estimating  the  target  position  T  by  that  value  rt  that 
maximises  the  probability  density  of  the  (ois,  u>R)-combinalion, 

For  any  event  of  combined  detected  positions  (ois  *  ti*R}  -  no  matter  how  far  apart  -  an 
estimated  target  position  t*  can  be  calculated  in  this  w  ay.  The  question  remains  how 
likely  it  is  that  both  contacts  result  from  the  same  target.  This  can  be  addressed  by 
calculating  the  total  probability  of  combinations  (Ws  ,WR)  that  are  equally  or  less  likely 
than  the  actual  observed  combination  (cos  ■  o>r>,  assuming  that  they  represent  the  same 
target  at  the  estimated  position  t*. 

Both  concepts  -  the  maximum  likelihood  estimate  of  the  target  position  and  the 
probability  of  the  observed  contact  combination  -  will  be  applied  first  in  a  one¬ 
dimensional  space  to  increase  our  understanding.  Thereafter  they  will  be  implemented 
for  two-dimensional  localisation  in  the  plane. 

2.2  One  dimensional  case 


2. 2,  /  Probability-density  Junctions 

Suppose  the  target  has  to  follow  a  very  narrow  sea  lane.  As  a  result  the  detection 
problem  becomes  one  dimensional,  making  the  analysis  rather  simple.  The  position  of 
the  contact  Ws  by  receiver  S  and  of  the  contact  W k  of  receiver  R  are  independent 
normally-distributed  stochastic  variables,  both  centred  at  the  target  position  T  but  each 
with  their  own  variances,  and  or  respectively. 
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That  i$T 

ws  ~  N(T,a*s) 

and 

wK  ~  N(T,<rl) 

The  probability-density  function  (pdf)  for  the  observed  combination  of  contacts 
(cos ,  (On)  becomes 


p{0)s,0)R)=  p{(Os)p{(Dfi)  = 


1  f  '  i 

- exp  — A 

lnosoR  \  2 


(3.1) 


with 


A  = 


a>s-T' 

- 

(  Vr-T  ) 

V  J 

1 

V  GR  ) 

(3.2) 


Remark  that  VV>  and  W&  are  stochastic  variables,  while  cuy  and  cur  represent  observed 
values. 


2,2,2  Maximum-likelihood  target  position 

For  a  given  realisation  (cos ,  cor)  the  value  of  p(a>s ,  ojh)  in  Equation  (3.1)  depends  on 
the  (unknown)  true  target  position  T.  Maximising p(cos, .  tu*?)  by  varying  T is  equivalent 
to  minimising  A ,  This  is  realised  if  7"  is  equal  to  r, .  the  ML  estimate  for  the  target 
position,  which  is  obtained  by  differentiating  Equation  (3.2)  and  is  given  by 


r.  = 


<r~s 


\ 


( 


\  °s  +aR  J 


K  + 


\a's  +(JR  J 


0)K 


(3.3) 


I,,  is  thus  a  linear  combi  nation  of  the  observed  contact  positions  cu.s  and  co#  and  is 
closest  to  the  observation  with  the  smallest  variance. 


2. 23  Probability  of  contact  association 

Derivation  of  the  probability  of  association  is  slightly  more  demanding  It  is  equal  to  the 
probability  of  combinations  of  contact  positions  {Ws ,  W#}  that  are  equally  or  less  likely 
than  the  actually  observed  combination  (cos  ,  cor).  under  the  assumption  that  the 
distributions  of  the  contacts  Ws  and  Wr  are  both  centered  at  i>.  In  the  following  we 
assume  that  ojr>  cos  ,  but  in  case  ojr<  ojs  a  similar  reasoning  can  be  followed. 


To  work  out  the  consequences  of  the  assumption  that  T  =  rt. ,  we  define  the  stochastic 
variable  as  follows 


T > 


K  + 


(3,4) 


r,..  the  ML  estimate  for  Tm  is  derived  from  the  contact  event  (cos*  wr)  and  therefore  is  the 
outcome  of  a  stochastic  process.  Taking  this  outcome  serious  imposes  a  restriction  on 
the  stochastic  contact  combinations  (  Ws ,  Wr)  taken  into  consideration  to  determine  the 
association  probability.  Only  those  combinations  are  allowed  that  yield  the  same  value 
rr  for  7\  The  assumption  is  thereby  transformed  into  the  condition  7V=  r*, ,  imposing  a 
relation  between  the  stochastic  variables  W$  and  . 
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With  this,  the  probability  of  contact  association  P„,„  becomes  equal  to  the  total 
probability  of  combinations  of  contact  positions  (Ws.W's)  that  are  equally  or  less  likely 
than  the  actually  observed  combination  ( co$ ,  o>r).  under  the  condition  that 
T'  =  re 


C  =  P{w, a a>„  A Ws  S ws  | T,  =  r, }  =  r-K  > a,  A A r,  r,} 


r,  =  r,  imposes  the  following  restriction  on  the  ( VVv,  WR)  combinations: 


r. 


<7's 


v  cr;  +  tx; 


(3.5) 


(3.6) 


Only  combinations  of  Ws  and  WR  that  satisfy  this  relation  can  result  in  the  estimated 
target  position.  Therefore  Ws  is  linked  to  via 


w  -(WMi 

s  ~~ 

oi 


(3.7) 


Consequently 

P{WK  >  coR  A  Ws  <  a 9S  A  Tt  =  r,}=  ^  W*  >  coR  a  IV,  =  (<T*  +  ~  a~xw*  |  (3.8) 


and 


P{T,=t,}=p\wi!Zr,AWs  = 


 (cr;  +  <7R)Tr  -  crsWH 
<JZ 


(3.9) 


Substituting  Equation  (3.7)  for  Ws  into  Formula  (3,1)  for  pi  WS,WR)  and  integrating 
over  WR  >  ojft  and  over  WR  >  xe ,  finally  gives  the  probability  of  association  of  the 
observed  combination  {a )$  ,ojr)\ 


p.,,<  =  I  -  erf 


a>R-a>s 

V2<7>-  +2(7  ■ 


E  J 


(3-10) 


where  erf(.)  denotes  the  error  function,  and  the  argument  of  the  error  function  is  equal 
to  the  so-called  Mahalanobis  distance  between  aj/fand  ws.  In  fact,  it  can  be  shown  that  if 
(Ws  ,  WR)  are  randomly  drawn  from  Ws  ~  N [Tr ,  <7* )  and  WR  ~  N(Tf>(JR)> 

respectively,  then  as  given  in  (3. 10)  is  equal  to  the  fraction  of  such  ( Ws  ,WR) 
combinations  that  are  more  distant  from  each  other  than  \a)s  - 

Figure  2. )  shows  the  effect  of  an  increasing  distance  between  the  observed  contacts  on 
their  association  probability.  We  will  give  some  numerical  examples  later  on  in 
Table  2.1. 


'  The  notation  is  given  in  ihe  list  of  symbols*  under  *  theory  general' 
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Figure  2. 1  Association 


2.3  Two-dimensional  case 


In  the  case  the  target  is  located  in  a  two-dimensional  plane,  the  target  position  T.  and 
the  contact  position  W*  and  WH  all  have  two  coordinates.  We  introduce  the  follow' ing 
notation 

Ws  =  {xs  y  ys )  W  R  =  {xR ,  y # )  and  T  =  ,  yT ) 

In  this  case  the  probability  distributions  of  the  two  detection  devices  become  bivariate 
normal  distributions,  both  centred  at  the  target  position  T,  but  with  their  ow  n  covariance 
matrices.  Cs  resp.  Cr  .  Le. 

Ws  -  yV(T,es)  and  VVK  -  N( TXR ) 

This  looks  rather  complicated,  but  it  is  not  if  it  is  presented  in  a  proper  way. 

Therefore  wre  first  present  multivariate  normal  distributions  and  then  go  back  to  the 
simpler  bivariate  case. 


23  J 


Intermezzo:  Multivariate  normal  distribution 

Suppose  v  is  an  n -dimens tonal  vector  w'ith  an  n- variate  normal  distribution. 
Then  the  pdf  of  v  =  (v/,  v?,  v„)T  is 


p(y)  = 


r expf“ ~ (v  -  M v  Y C (v-Mf) 


V(2^)"Dff(  C()  ^  2 

w  ith  p*  the  vector  of  expected  values  of  v,  CT  the  covariance  matrix  of  v 
<*ll  ^12  -  CT]„1 

CT2l  <J22  &in 


C  = 


_&tt\  &fl2  &nn  J 


(3.11) 
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with 

<Ttm  =  ElV*  "A,  )T  (Vm  "  Mvm  )}=  CFmk 

and  Det(Cy)  the  determinant  of  Cv  * 

We  define 

-  yj&kk 
and 


a. 

A*  =  ~rzr 

A  A  I 


For  the  case  n  =  l  the  well-known  pdf  of  the  one-dimensional  normal  distribution 
results. 


23.2  Probability-density  fimctions 

The  positions  \VS  detected  by  receiver  S  have  the  following  bivariate  normal 
distribution  (/i  =  2): 

U.s-  N{{xT,yT)Xs) 


with 


Cs  = 


(  _2 

ox 

a12  ' 

f 

— 

Kai\ 

V 

*XS 

^XS  ®YS 


Ps  &  X5  &YS 


YS 


\ 

) 


There  is  a  change  in  notation  here:  the  dimensions  1  and  2  are  relabelled  by  X  and  Y 
respectively. 

Inserting  the  covariance  matrix  Cs  into  the  pdf- Formula  (3.11)  results  in 


p{v>s)=  p{xs.ys)  = 


2ko  ys  \\  Ps 


exp 


— -A* 


(3.12) 


with 


1 

'xs -xT  ' 

1 

-  2Ps 

(xx  -  XT  ) 

( ys  ~yT ' 

+ 

fys-yr  V 

1 

V  °*JfS  J 

v  &XS  J 

{  <*YS  J 

v  J  J 

(3.13) 


For  the  positions  VV'h  observed  by  receiver  R  each  subscript  'S’  must  be  replaced  by  an 
‘R\  Since  the  stochastic  contact  positions  VV.S  and  WK  are  assumed  to  be  independent, 
the  pdf  for  a  combination  of  both  contacts  is  the  product  of  their  individual  pdf  s. 

For  the  measurement  (tos ,  cur)  =  ((xs  ,ys ) .  (Xr  .>’«))  the  probability  density  becomes 


pta.ys.  *«.>’*)= 


xs  &ys  v/l  ~Ps  */i  ~Pr 


1  1 - -expf~l(Av  +  (3, 

cry„Jl-p;  \  2  ) 


14) 


Remark  that  the  assumption  that  both  contacts  come  from  the  same  target  is  already 
incorporated  in  the  ,4  =  4s  +  ,4*  pan  of  this  formula.  The  target  position  (Xr.yr)- 
however,  is  yet  unknown  and  must  be  estimated. 
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233  Maximum  -likelihood  ta rget  position 

For  a  given  combination  (tuHs  ,  gir)  the  value  ofp(ti)s  *  cur)  depends  on  T  =  (xy-,  Vy). 
Maximising  p{ tus  »  u>r)  by  varying  T  is  equal  to  minimising  A.  Making  both  partial 
derivatives  of  A  with  respect  to  Xj  and  yy  equal  to  zero,  results  in  the  ML  estimate  for 
the  target  position  te  =  (r*  .  ry),  wrhich  can  be  shown  to  be  equal  to 

tc  -  cos  +  M  (to  R  -  (0S )  (3, 1 5) 


with 


M  = 


r) 

P) 


and 


a  = 


_  ®xs&3ps  0  ®  XS^YR  PuPs^f  XS&YS® XK^YX 


B 


(3.16) 


P  = 


_  G'xs&rsips  0  &XF^YS  PrPs^XS^YS^XR^YR 


B 


7  = 


2 

Pr^XR^YR^JCS  ~  P xsO ysff XR 

B 


3  = 


Pr^  XR^YR^YS  PsQjsQjxOjR 

B 


B  —  O xs(7ys{fis  l  )  +  ^XR^YR^Pr  0 

1  2  2  2  *y 

—<JxS  ®YR  ~  &XR&YS  +  ~ Ps  Pr &XS  &YS &XR &YR 

Usually  tv  will  not  just  be  a  linear  combination  of  the  observed  contact  positions  cus  and 
(Ur,  because  M  ^  «J  .  If  a  ~  fi  and  y  -  S  =  0.  (ix  ,  ry)  will  be  on  the  line  between 
(xs  -  ys)  and  (xj?  *  vr)*  In  this  case 

=  and  ps  =  pR 

&YS  &YR 

and  the  pdf-contours  of  Ws  and  VVR  will  have  exactly  the  same  shape  and  orientation* 
but  can  have  a  different  scaling  factor  This  is  demonstrated  in  the  upper  panels  of 
Figure  2,2,  where  the  target  position  is  on  the  line  between  the  two  observed  contacts. 
The  fact  that  the  ML  estimate  of  the  target  position  te  is  usually  not  on  the  line  between 
the  two  observed  contact  positions  tu.s  and  u)r*  is  an  unexpected  result  that  is  important 
for  target  localisation.  This  is  illustrated  in  the  lower  panels  of 
Figure  2*2,  By  plotting  the  probability  density  contours  of  the  contacts*  these  panels 
give  an  impression  of  how  the  probability  distributions  of  Wg  and  Wr  can  influence  the 
ML  target  position.  In  Table  2. 1  the  input  and  output  of  Figure  2,2  is  presented 
numerically,  Note  that,  in  each  panel*  the  contours  shown  for  and  cor  do  in  general 
not  correspond  to  the  same  pdf- value*  hut  are  the  ones  that  pass  through 
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Table  2. 1  Numerical  inpuls  and  outputs  used  to  generate  the  four  panels  of  1'igure  2.2  (arbitrary  units). 


upper  left 

upper  right 

lower  left 

tower  right 

input 

*s 

2000 

2000 

2000 

2000 

ys 

3000 

3000 

3000 

3000 

Oxs 

15 

15 

20 

20 

Gys 

20 

20 

40 

40 

Ps 

0.4 

O.S 

0.4 

0,8 

Xr 

2030 

2030 

2050 

2050 

Vfl 

2960 

2960 

3000 

3000 

OXR 

30 

30 

40 

40 

Oyr 

40 

40 

60 

60 

PR 

0.4 

0,8 

-0.4 

-0.8 

output 

Xr 

2006 

2006 

2011.4 

2016.5 

yr 

2992 

2992 

3013,5 

3030.7 

association 

probability 

0.24 

0.02 

0.59 

0.60 

0 

0.2 

0.2 

0.23 

0.33 

p 

0.2 

0.2 

0.34 

0,46 

Y 

0 

0 

0.09 

0.20 

6 

0 

0 

0.27 

0.61 

fe4{MU*|.  •glrwft  iJMKR) 


i»0  1BS0  2000  2050  2100  2150 


Sjirrik 


194C 


IBM  1  too  2000  2020  20*0  MM  30M  HOC  2120 


TB0Q  ISfiD  ?ODO  7060  JIM 


•R  |WW#I.  (W»CHS 


Figure  2,2  Example  pdf  contours  of  contact  positions  W*  and  Wr  centred  at  Ms  (red  doi )  and  iur  ( blue 
dot),  wiih  maximum  likelihood  target  position  Te  (black  dot)-  Upper  panels:  target  tc  on  line 
between  m  and  cor  :  lower  panels:  targei  tv  not  on  the  line  between  Us  and  or  -  Parameter 
values  given  in  Table  2.1, 
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2.3, 4  Probability  of  contact  association 

As  in  the  one-dimensional  case,  the  probability  of  association  is  set  equal  to  the 
likelihood  of  getting  the  observed  contact  positions  cos  and  tuR  (or  less  likely  values) 
from  the  distributions  of  Ws  and  WR  ,  under  the  assumption  that  these  distributions  are 
both  centred  at  t*.  Let  us  define  the  stochastic  variable  T*  as 

Te  =  Ws  +  M  ■  ( WR  -  Ws )  (3.17) 

Following  the  same  reasoning  as  in  the  one-dimensional  ease,  the  assumption  that 
T  -  te  is  transformed  into  the  condition  T*  =  t*  ,  which  imposes  the  restriction  on  the 
stochastic  contact  combination  (Ws,Wr)  that  each  combination  should  yield  the  same 
ML  estimate  t*  for  T  as  the  observed  combination  (ois  .  gjr)  does. 

With  this  the  probability  of  contact  association  becomes  equal  to  the  total  probability  of 
combinations  of  contact  positions  (Ws  ,  WR)  w  hich  are  equally  or  less  likely  than  the 
actually  observed  combination  (ms  *  cdr)  under  the  condition  that  Te  =  t*  . 

In  the  tw  o-dimensional  case,  however,  the  definition  of  Mess  likely"  has  to  be  adjusted. 
Here  we  call  positions  Ws  less  likely  than  the  observed  value  ais 
if  pdfs(Ws)  S  pdfs(u>sX  with  pdfs(iDs)  the  probability  density  at  ids-  So  less  likely 
positions  Wg  are  all  points  outside  the  probability  density  contour  going  through  w s* 
which  can  be  shown  to  be  an  ellipse.  We  denote  all  vectors  inside  this  ellipse  by  £(cds). 
The  same  can  be  done  tor  the  R  system,  by  replacing  all  subscripts  'S'  by  hR\ 

With  all  this  background  the  association  probability  PiiSS  is  defined  as 
PISS  =  P{ WK  i  £(wk  )  V  Ws  e  E{ ros  )|T(  =t,}= 

/>{\Vrg  £(o)r)aWs€  £(cos)ATe  =tJ  (3.18) 

H t.=T.} 

The  condition  Te  =  x*  is  equivalent  to  the  following  expression  for  W* 

Ws  =  (I  -  M rl  ■  (t.  -  M  ■  WR  )  (3.19) 

With  this  we  get 

p{ WR  e  E{ to JaWsg£((0s}aT(=tJS 

p{\\r  <=E{ o>R )  a  Ws  €  E{ cos )  a  Ws  =  (I  - M )-’  ■  (t.  - M  ■  VVR )} 

and 

P{T  =  tJ=p{\Vr  g  R2  a  Ws  =  (i  - M)~*  •  (t.  - M  •  WR )} 

wfhere  R:  denotes  the  entire  two-dimensional  plane. 

Given  the  value  of  xtf  =  (r*  ,  r^)  we  redefine  (xr  ,  Vr)  =  (Xx  »  Tk).  Making  (Xr  ,  »)  the 
integration  variable  running  over  £(cdr)  and  setting 

U, ,  ys )  =  (I  -  M r1  {{xT ,  yT  )r  -  M{xh  ,  yft  f )  (3.20) 
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Equations  (3. 12  -  3. 14)  can  be  used  to  get 
P{ WR  e  E{(o „  ) a  Ws  e  E( o»s )aT,=tJ= 


J  J  MrMs  expf  -  i  (xR ,  )  +  As  (_xs ,  v5  ))  ps dys  ±\B  dy K  = 


MrHs  |  eXp[-^R(jCR-.VR)  {  expf- ^  AS  (*5  ’  >S  )  Vw.Vj 


(3.20) 


where 

Mtt= - 1  f—j  0-21) 

XR&YR  \  ^  “  Pr 

and  a  similar  expression  for  jUs  <  and  where  As  and  AR  are  defined  in  Section  2,3,2. 
Remark  that  (x$  ,  y$)  depends  on  (xr  ,  \r).  Although  the  fact  that  Cr  and  Cs  are 
symmetric  matrices  allows  for  several  nice  tricks  (such  as  applying  basis 
transformations),  the  second  integral  results  in  a  complicated  error-function  expression. 
Instead  of  setting  out  on  the  heroic  task  of  analytically  integrating  this  error  function 
(certainly  prize  wanning  if  successful),  we  have  chosen  to  determine  this  probability  in  a 
numerical  way. 

23  A,  I  Alte  nut  ti  ve  defirt  i ti o n  of  pro  hah  i l i rv  of  c on  tac  t  association 

Remark  that  an  alternative  definition  of  the  association  probability  in  Equation  (3. 18) 
would  be  possible.  It  could  also  have  been  defined  by 

p{ W„  C  £(©  B )  A  Ws  g  E{ms  )  \  X  =  Tt } 

what  leads  to  an  equal  or  smaller  association  probability  .  This  choice  cannot  be  made  in 
one-dimensional  situations,  where  a  less  likely  value  Wk  >  Mr  can  only  be  combined 
with  a  value  Ws  <  cos  (also  less  likely  than  the  actual  observed  value  ms)  if  the  condition 
Tr  -  xe  must  be  met.  For  contacts  in  a  2-D  plane,  however,  it  is  possible  that  a  less  likely 
R -contact  position  (Wr  outside  the  ellipse  E{ cor))  has  an  associated  value 

ws=(i-m)-'(t,-mwb) 

that  is  inside  the  ellipse  E(cos)  and  is  therefore  more  likely  than  e)s.  Here  we  have 
chosen  to  define  a  contact  combination  (Wr  ,  Ws)  less  likely  if  at  least  one  of  those 
values  is  outside  its  pdf  ellipse.  This  choice  gives  the  highest  association  probability. 


ONGERUBRICEERD 


ONGERUBRICEERD  |  TNO  report  |  TNO-DV  2006  A518 


21/51 


3  Understanding  the  association  probability 


All  the  relevant  Equations  presented  in  Chapter  2  are  implemented  in  Mat  lab  to  solve 
them  numerically.  It  is  worthwhile  to  analyse  some  outcomes  to  understand  the 
processes  that  determine  the  association  probability.  First  the  successive  steps  for 
calculating  the  association  probability  are  illustrated  in  Figure  3.1. 

Thereafter  the  results  of  a  carefully  chosen  set  of  contact  situations  are  presented  and 
compared  with  the  one-dimensional  solutions. 

3.1  Visual  presentation  of  the  association  probability 

For  the  situation  given  in  the  lower  right-hand  panel  of  Figure  2.2,  the  upper  left-hand 
panel  of  Figure  3. 1  shows  the  contacts  Or  =  (2050 , 3000),  tils  =  (2000 . 3000)  with  the 
estimated  target  position  r,=  (2016.5 . 3030.7)  and  the  ellipses  £(or)  and  £(os). 

For  all  Wr  values  inside  £(Wr)  with  ,t«  =  1990  the  coupled  W.s  values  are  calculated. 
Some  of  these  Ws  values,  however,  are  outside  their  £(Os)  ellipse. 

In  the  upper  left-hand  panel  only  those  27  combinations  of  Wr  and  Ws  values  are 
shown  that  are  both  inside  their  respective  ellipses.  For  this  reason  some  Wr  and  Ws 
values  seem  to  be  missing. 

Each  Wr  has  a  pdf  value  according  to  its  bivariate  normal  distribution.  Its  coupled  Ws 
value  has  its  own  pdf  value.  The  logarithm  of  these  pdf  values  and  the  logarithm  of  their 
product  are  shown  in  the  upper  right-hand  panel  of  Figure  3. 1  for  each  of  the  27 
(Wr  .  Ws)  combinations. 

In  the  lower  left-hand  panel  of  Figure  3.1  for  each  possible  Wr  position  in  the  plane 
(whether  inside  or  outside  £(ojr))  the  coupled  Ws  value  is  determined  and  the  product 
of  their  pdf  values  is  assigned  to  the  W'r  position.  This  shows  the  conditional 
probability  densities  for  this  specific  case. 

In  the  lower  right-hand  panel  of  Figure  3.1  all  probability-density  values  for  which  the 
(Wr  .  Ws)  combination  is  inside  £(u)r).  respectively  £(a>s).  are  set  equal  to  0.  The  total 
probability  density  present  in  this  panel  (viz.  3.5e-5)  divided  by  the  total  probability 
density  present  in  the  lower  left-hand  panel  (viz.  5.8e-5)  gives  a  probability  of 
association  of  0.60  for  the  contacts  Or  =  (2050 . 3000)  and 

o>s  =  (2000 . 3000)  (for  the  numerical  input  values  used:  tixs  =  20.  rryS  =  40.  ps  =  0.8. 

-  40.  frVR  =  60.  pR  =  -0.8). 
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Figure  3  J  Upper  left:  the  dots  show  the  observations  wr  (blue)  and  ta*  (red)  and  the  most-likely  target 

position  tc  (black).  The  arrows  show  the  coupling  of  Wr  and  Ws  positions  inside  the  ellipses 
E(g>r)  respectively  £(tus).  Upper  right:  probability  density  values  of  (Wr  t  Ws)  combinations 
around  xr  Middle  and  lower  panels:  colour  maps  of  pdf-values  contributing  to  probability  of 
association:  the  probability  equals  Lhe  sum  oF  the  pdf- values  in  the  right-band  panel  divided  by 
the  same  sum  of  the  left-hand  panel.  Hie  middle  panels  are  a  zoom  of  the  lower  panels,  at  the 
same  scale  as  the  upper  left-hand  panel. 

3.2  Practical  way  to  calculate  the  association  probability 

The  above  illustration  suggests  the  following  steps  for  lhe  calculation  of  the  association 

probability. 

1  For  a  pair  of  observed  contacts  (cor  ,  cos)  determine  the  ML  target  position  t*.  using 
Equation  {3. 15). 

2  Calculate  the  ellipses  E(cor)  and  £(cos)  from  toK.  i&s  and  t*. 

3  Choose  a  sufficiently  large  integration  area  around  t *  in  the  two-dimensional  plane. 
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4  For  each  point  \VK  in  this  integration  area  determine  the  accompanying  VVS  position 
that  combines  with  Wr  to  T*  =  t*. 

5  For  each  (\VR  ,  Ws)  combination  calculate  (he  probability-density  values  pdf(WR) 
and  pdf(Ws)t  using  Equation  (3.12)  and  their  product  p(WR)  =  pdf(WRf  pdf(  VVS }. 

In  this  way  a  'probability  density1  value  p( Wr)  is  assigned  to  each  point  Wr  in  the 
integration  area. 

6  Integrate  p(WR)  over  the  integration  area. 

7  Integrate  p( \VR)  over  that  part  of  E( toR)  for  which  the  accompanying  Ws  values  are 
also  inside  E(gj s)- 

8  The  association  probability  for  die  observed  contact  combination  (iur  ,  tus)  is 
1  minus  (outcome  of  7)  divided  by  (outcome  of  6). 

3,3  Analysis  of  association  probability 

To  study  the  mechanisms  of  the  association  process,  a  set  of  mutually  comparable 
situations  is  considered.  The  contact  positions  are  chosen  as  follows:  In  all  cases 
x$  =  2000,  \s  =  3000  and  }r  =  3000.  Only  the  xr  coordinate  is  varied,  from 
=  *s  +  &xrI 4  (via  xr  -  Xs  +  <tXr/2  ,  Xr  —  Xs  +  3rr^/4  ,  xr  =  xs  +  )  to 

Xr  —  Xs  +  2(7xr 

For  the  distributions  of  the  contacts  the  following  choices  are  made: 

1  qXr  =  2oxs 

2  oyr = 2oys 

3  either: 

a  pR  =  ps  =  0.8  (so  y  =  5  =  0  and  Te  on  line  w  ith  ior  and  u>s) 

b  Pr  =  -p$  and  ps  =  0.8  (so  y8  *  Q  and  ie  is  off  line  with  o>H  and  ©s) 

The  values  of  a xs  and  aYs  are  varied  (see  Table  3, 1 }.  The  first  two  choices  mean  that 
a  ~  /(,  see  Section  2.3.3. 

Table  3,1  summarises  the  results  of  this  investigation.  From  this  Table  the  following 
conclusions  can  be  drawn: 

1  The  three  cases  (oXs  =  20,  ays  =  40),  (fJXs  =  200,  ays  -  4)  and  {a Xs  =  2,  ays  =  400) 
only  differ  by  a  scaling  factor  in  both  coordinate  directions,  but  all  calculated 
probabilities  are  exactly  the  same.  Therefore  it  suffices  to  limit  our  attention  to  the 
first  case  ( aXs  =  20,  ays  -  40). 

2  The  association  probabilities  for  differently  oriented  contact  distributions 
(that  is  when  pa  -  -ps  and  so  t*  is  not  on  the  line  between  and  tps)  are  much 
higher  than  for  distributions  w  ith  the  same  orientation. 

3  For  distributions  w  ith  a  different  orientation  (pR  =  -ps)  the  association  probability  is 
considerably  higher  than  could  be  expected  on  the  basis  of  a  one  dimensional 
analysis.  But  if  both  contact  distributions  have  the  same  orientation  {j)r  -  ps)  then 
the  association  probability  is  higher  than  the  one-dimensional  probability  for  two 
near  contacts,  but  lower  for  twro  distant  contacts.  Further,  ps  =  0,8  is  a  high  value. 
For  lower  values  of ps  the  association  probabilities  in  case  both  distributions  have 
the  same  orientation  {jjr  =  ps)  increase,  but  remain  lower  than  for  differently 
oriented  distributions  (/)r  =  «/>$). 

To  conclude,  the  more  the  distribution  of  the  S  contact  differs  from  the  distribution  of 
the  R  contact,  the  more  important  the  contribution  of  the  two-dimensional  analysis 
becomes. 
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Table  3. 1  Results  of  associating  i wo  contacts  with  different  parameters.  See  text  above  Table  for  details. 


xc  on  line  between  wk  anti  (pfl  = 

Ps) 

tc  noi  on  line  between  ojk  and  ojs  (/>k  =  -Ps) 

x*  -  = 

<jyJ4 

3<W4 

aXK 

<W4 

<W2 

3<W4 

°\ R 

_ 

°XS  —  20,  Oys  -  40 

** 

2010 

2020 

2030 

2040 

2080 

2010 

2020 

2030 

2040 

2080 

XT 

2CKJ2 

2fMJ4 

2006 

2008 

20!  6 

21X13.6 

21X17.2 

2010.8 

2014.4 

2020 

Vj 

3000 

3000 

3000 

3000 

3000 

3006.7 

3013.3 

3020 

3026,6 

3053,2 

ass  prob 

0.93 

0.76 

0.54 

0.33 

0.01 

0.98 

0.91 

0.82 

0.69 

0.24 

1-dim  prob 

0.82 

0.65 

0.50 

0.37 

0.07 

0.82 

0.65 

0.50 

0.37 

0,07 

(fxs  =  200,  an  =  4 

Xft 

2100 

2200 

2300 

2400 

2800 

2100 

2200 

2300 

2400 

2800 

XT 

2020 

2040 

2060 

2080 

2160 

2036 

2072 

2108 

2144 

2288 

Vr 

3000 

3000 

3000 

3000 

3000 

3000.7 

3001.3 

3002 

3002.7 

3005.3 

ass  prob 

0.93 

0.76 

0.54 

0.33 

0.01 

0,98 

0,91 

0.82 

0.69 

0.24 

bdim  prob 

0.82 

0.65 

0.50 

0.37 

0.07 

0,82 

0,65 

0,50 

0,37 

0,07 

oxs  =  2,  (7^  =  400 

2001 

2002 

2003 

2004 

2008 

2001 

2002 

2003 

2004 

2008 

-*r 

21KKJ.2 

2(HX).4 

21XX1.6 

2000,8 

2001.6 

2000.4 

2000*7 

21X11.1 

21X11. 4 

21X12.9 

Vr 

3000 

3000 

3000 

3000 

3000 

3066,5 

3133 

3199.5 

3266 

3532 

ass  pmb 

0.93 

0,76 

0.54 

0.33 

0.01 

0.98 

0.91 

0.82 

0.69 

0.24 

1-dim  prob 

0.82 

0.65 

0.50 

0.37 

0.07 

0.82 

0.65 

0,50 

0,37 

0.07 
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4  Performance  estimation  by  means  of  simulations 


41  Introduction 


in  order  to  study  the  consequences  of  the  theory  exposed  in  Section  23,  and  to  evaluate 
different  data-association  rules,  we  have  simulated  two  data  sets,  corresponding  to  the 
outputs  of  the  two  receiver  systems  S  and  R.  The  detection  performance  of  S  and  R 
separately  can  then  be  compared  to  that  of  the  combined  system.  This  comparison  is 
done  using  classical  ROC  curves. 

Section  4,2  describes  the  simulations  that  are  used  for  this  analysis.  In  Section  4,2.3  the 
data-fusion  algorithm  is  described  in  detail.  Section  4.3  presents  the  simulated  scenario. 
In  Section  4.4  the  results  for  the  single  sonars  will  be  given.  Sections  4.5  to  4. 10  present 
die  results  after  combining  or  fusing  both  data  sets. 


Before  describing  the  simulations  and  their  results,  it  is  important  to  explain  a  few 


terms: 


‘detection*  or  ‘contact*: 
‘target*: 

'target  detection’: 

‘false  alarm*: 

‘combinina'  data  sets: 


General  threshold  crossing. 

‘real’  target  inserted  in  the  simulations. 

Threshold  crossing  corresponding  to  a  Teal'  target. 
Threshold  crossing  not  corresponding  to  a  real  target*  but 
due  to  a  statistical  fluctuation  of  the  background. 

The  contacts  in  the  (two)  original  data  sets  are  gathered  in 
one  list. 


*  'fusing'  data  sets:  The  contacts  in  the  (two)  original  data  sets  are  gathered 

by  applying  a  fusion  algorithm  to  them  and  then  making  a 
list  containing  all  contacts  after  this  fusion, 

‘Combining’  the  data  sets  can  thus  be  seen  as  'fusing'  them  without  any  fusions  to 
occur. 


42  Simulations 


4.2. 1  Da  ta  -s  im  it  la  tion  proce  cl tt  re 

In  the  simulations  we  have  assumed  two  sonars  to  be  present,  a  monostatic  one 
(sonar  i )  and  a  bistatic  receiver  (sonar  2),  Both  sonars  are  characterised  by  their 
position,  and  by  a  set  of  uncertainties  in  the  follow  ing  relevant  quantities:  sound  speed, 
position,  receiver  heading,  bearing  estimates,  time  stamps. 

These  uncertainties  represent  the  errors  that  are  likely  to  be  present  in  real  sonar 
systems,  as  described  e.g.  in  [I]. 

The  simulations,  used  to  assess  the  benefits  of  the  data-association  algorithm  presented 
in  the  previous  Chapters,  tire  set  up  as  follows: 

1  General: 

a  Define  the  uncertainty  in  the  sound  speed  (<Jc),  This  value  is  assumed  to  be 
equal  for  both  sonars  and  constant  over  pings,  and  can  be  estimated  from 
experience  with  sea- trial  experiments. 

b  Two-dimensional  propagation  is  assumed,  i.e.,  no  depth  information  is 
included. 
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2  Sonars: 

a  For  sonar  i  {i-  L  2)  define  its  position  (xV  ( ,  yu). 

b  For  each  sonar,  define  the  sampled  beams. 

c  Define  the  sample  frequency  (assumed  equal  for  both  sonars)  and  the  ping 
repetition  rate.  The  latter  defines  the  covered  area  for  both  sonars. 

d  For  each  sonar,  generate  sonar  signals  for  each  (range  .  bearing)  cell. 

The  sonar  signals  are  intended  to  represent  real  sonar  signals  after 
beamforming  and  matched  filtering,  as  in  practice  detection  is  applied  after 
these  two  processing  steps,  The  signals  in  the  simulations  are  drawn 
randomly  from  a  Rayleigh  distribution,  as  this  is  the  distribution  that  is 
most-often  assumed  to  characterize  real  data. 

e  For  sonar  i  (i  =  ! ,  2),  define  the  errors  in: 
i  Position  (GxXJ  ,  GyStt ,  GxySti). 

h  Receiver  heading  (c#). 

iii  Bearing  (G^), 

iv  Time  synchronisation  (GU« 

3  Targets: 

a  Define  a  number  of  targets  nt  at  fixed  positions  (x{J  f  y^),  j  =  1 . nt. 

b  Define  the  target  signal -to- noise  ratios  SNRjj  of  target  j  with  respect  to 

sonar  L  These  signal -to- noise  ratios  relative  to  both  sonars  can  be  different 
due  to  the  fact  that  both  sonars  see  the  target  at  a  different  aspect  angle,  for 
example. 

c  Determine  the  theoretical  uncertainties  (nx$j  *  Gvfj  .  a xy,j)  in  the  target 
positions,  based  on  (Gc  ,  gxsj ,  0\\j  ,  GxysJ  ,  .  o<pt ,  Off),  For  this,  the 

Equations  of  [3]  are  used.  Note  that  the  expected  uncertainty  in  the  target 
position  will  in  general  be  different  for  sonars  1  and  2,  as  the  parameters 
(G.xs  t ,  avy, ,  GxySri ,  ,  ott)  themselves  can  be  different  for  both 

sonars.  In  addition,  the  target  will  in  general  be  at  a  different  relative 
position  relative  to  the  sonars. 

d  For  sonar  /  and  for  each  ping  p,  generate  an  ’observed'  target  position 
(Xtjxp  *  based  on  the  true  target  position  (.t,r; ,  ytJ)  and  (Gxui ,  avM  , 
Gxyu). 

4  Data: 

a  For  both  sonars  and  for  each  target,  find  the  (range,  bearing)  cell  Cup  that 
is  closest  to  the  'observed'  target  position  (j itj^p  ,  ytjxP). 

b  At  the  cells  C(jtPf  generate  new  sonar  signals  with  signal-to-noise  ratio 
SNRjj, 

c  Take  the  logarithm  of  the  data  in  every  (range,  bearing)  cell. 

4,2.2  5 ing  le -sonar  de  tec  lion  p  rocedu  re 

After  generating  the  data,  as  described  in  the  previous  Section,  the  detection 
performance  of  both  single  sonars  is  evaluated  in  the  following  way: 

For  both  data  sets  /,  apply  amplitude  thresholds  At  running  from  some  low  number 
(e.g.  9  dR)  to  some  high  number  (e.g.  30  dB). 

For  each  threshold  At,  count  the  number  of  threshold  crossings  in  all  cells  CfJip 
This  number  represents  the  number  of  target  detections  NiUtP  (Ar)  for  sonar  i. 

For  each  amplitude  threshold  At*  count  the  total  number  of  threshold  crossings  outside 
cells  This  number  represents  the  number  of  false  alarms.  N^p  (Ar)  for  sonar  f. 

The  probabilities  of  detection  and  false  alarm  are  determined  by  measuring,  for  each 
ping  p,  the  number  of  detections  Nsi  iiP(AT)  and  false  alarms  Nf(lxp(Ar)  as  described 
above.  These  numbers  are  then  summed  over  ail  pings,  which  gives  the  total  number  of 
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detections  and  false  alarms  over  all  pings.  These  total  numbers  are  then  divided  by  their 
respective  number  of  opportunities  to  have  a  threshold  crossing.  For  the  number  of 
detections*  Nt{  iif}(ATh  the  number  of  opportunities  is  equal  to  the  number  of  targets  per 
ping  times  the  number  of  pings.  For  the  number  of  false  alarms,  N/^piAr)*  the  number 
of  oppon unities  is  equal  to  the  number  of  (range  *  bearing J-ce  11s  minus  the  number  of 
targets*  times  the  number  of  pings. 

In  summary*  the  probabilities  of  detection  and  false  alarm  for  sonar  i  (i  -  l*  2)  are 
determined  via 


P*AAr)  = 


N  . 

pings 


l 

#  targets 


p= i 


(A,.) 


(4. 1 ) 


P  fa. i  (Ar  )  - 


V  .  # (range, bearing)cells  -  #targets 


I 


p=i 


(4.2) 


Plotting  pd.  {At  )  versus  pfa  i  (Af  )  generates  a  so-called  ROC  curve  [6]*  This  is  a 
well-known  and  well-accepted  measure  for  estimating  the  sonar  performance. 


4. 2, 3  Data- fusion  procedure 

After  having  identified  the  two  sets  of  threshold  crossings*  for  every  amplitude 
threshold*  these  two  sets  are  fed  to  the  fusion  process.  In  practice,  this  fusion  algorithm 
is  implemented  as  follows: 

1  The  input  data  sets  Ci  and  C2  consist  of  all  contacts,  above  the  relevant  threshold* 
observed  by  sonar  1  and  2?  respectively*  Each  contact  in  the  data  sets  has  seven 
labels: 

a  (x  *  y).  denoting  the  contact  position  (in  metres)* 

b  (o.v .  ay  *  <xry),  indicating  the  expected  uncertainty  in  the  contact  position 
(see  Section  4.2.1), 

c  SNR,  containing  the  mean  contact  signal -to-noise  ratio, 
d  ID,  a  label  that  specifies  if  the  contact  is  a  false  alarm  {ID  =  0)*  or 
corresponds  to  any  of  the  targets  {ID  =  1,  .  nt). 

To  simplify  things  in  the  algorithm  implementation*  Q  is  assumed  to 
contain  not  more  elements  than  C2.  If  it  does,  C]  and  C 2  are  interchanged. 

2  Calculate  the  association  probability  between  any  pair  of  contacts  from  both  data 
sets.  The  algorithms  developed  in  Section  2.3  are  used  for  this. 

3  For  each  contact  q  in  data  set  Q  *  find  the  contact  q  in  data  set  C2  for  which  the 
association  probability  Pass  is  largest. 

4  Check  association  by  drawing  a  random  number  between  0  and  L  If  this  number 
is  smaller  than  Pa„*  the  two  contacts  cj  and  cy  are  fused 
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5  Each  contact  that  exists  after  the  fusion  process  is  passed  to  the  output  data  set 
and  has  the  following  tour  labels: 

a  (x/f  \f),  the  ML  contact  position  (in  metres),  as  derived  in  Section  2.3. 
b  SNRj ,  the  mean  contact  SNR,  For  fused  contacts  it  is  the  maximum  of  the 
two  SNRs  of  the  original  contacts.  For  non-fused  contacts  it  contains  the 
SNR  of  the  original  contact. 

c  IDf,  a  label  that  specifies  the  status  and  origin  of  the  contact.  The  status 
can  be: 

i  ‘not  fused' 

ii  ‘fused’ 

6  For  non-fused  contacts  the  contact  can  originate  from 

i  a  true  target 

ii  a  false  alarm 

7  For  fused  contacts  the  contact  can  originate  from: 

i  twice  the  same  true  target  (e.g.,  target  2  as  observed  with 
sonar  1  and  target  2  as  observed  with  sonar  2), 

ii  twfo  different  true  targets  (e.g.  target  3  as  observed  with 
sonar  1  and  target  5  as  observed  with  sonar  2). 

iii  a  true  target  and  a  false  alarm. 

iv  two  false  alarms. 

It  is  necessary  to  keep  track  of  all  these  possibilities  to  make  sure  that  the  data- 
fusion  process  does  connect  only  valid  pairs  of  contacts  and  not  spurious  close 
pairs. 

There  is  one  free  parameter  in  the  fusion  process  as  we  have  implemented  it:  the  size 
and  shape  of  the  covariance  ellipses  around  the  contact  positions  are  used 
(see  Section  2.3.2)  to  calculate  the  association  probabilities.  In  principle,  these  ellipses 
should  be  the  same  as  those  that  are  used  to  generate  the  ‘observed'  target  positions 
(see  Section  4.2.1).  However,  wre  have  allowed  the  ellipses  used  for  the  fusion  process 
to  be  F times  larger.  The  larger  F^Up^  .  the  easier  targets  are  associated. 

The  introduction  of  this  free  parameter  seems  somewhat  ad  hoc.  but  will  turn  out 
to  be  useful  later  on. 

4.2.4  De  tec  ti on  procedu  re  after  comb  in  ing  da  ta  se  ts 

After  combining  the  two  original  sets  of  contacts,  one  data  set  remains  that  contains  all 
contacts  of  the  two  original  data  sets.  In  this  data  set: 

•  True  detections  are  identified  as  the  true  detections  of  the  original  data  sets, 

*  False  alarms  are  identified  as  the  false  alarms  of  the  original  data  sets. 

It  should  be  ensured  that,  if  a  target  is  observed  by  both  receivers,  it  is  not  counted 
twice  as  a  true  detection.  That  is,  per  ping  the  number  of  true  detections  can  never 
exceed  the  number  of  targets  present  in  the  simulations. 

This  wray  of  combining  the  data  sets  w  ill  be  referred  to  as  ‘option  1 '  hereafter, 

4.2.5  Detection  procedure  after  fusing  data  sets 

After  fusing  the  twro  original  sets  of  contacts,  one  data  set  remains.  For  this  data  set,  we 

have  defined  different  options  which  each  have  a  different  ‘definition1  of  true  detections 

and  false  alarms.  These  definitions  are  summarised  in  Table  4.  L 

For  example,  in  option  3  a  contact  after  fusion  is  labelled  a  true  target  if  it  consists  of 

either: 
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•  A  fused  contact  that  originates  from  two  identical  true  targets  (T,T,). 

•  A  (non -fused)  contact  that  consists  of  a  true  target  (T,). 

•  A  fused  contact  that  originates  from  a  true  target  and  a  false  alarm  (T|  F). 

•  A  fused  contact  that  originates  from  two  non-identical  true  targets  (T,Tt). 

Options  2  and  3  can  be  seen  as  variants  of  the  ‘OR'  rule:  if  one  of  the  receivers  detects 
the  target  it  is  taken  into  account.  Options  4  and  5  can  be  seen  as  variants  the  'AND' 
rule:  both  receivers  have  to  detect  the  target,  and  these  two  detections  should  fuse  with 
each  other. 

For  the  labelled  data,  and  for  each  amplitude  threshold  At.  the  numbers  of  true  targets 
and  false  alarms  with  SNR  >  Ar  are  then  counted.  For  options  12  -  15  the  numbers  of 
true  targets  and  false  alarm  with  SNR  >  Ar  +  AAr  are  counted. 

Here  AA?-  is  the  increase  in  threshold  that  is  applied,  which  is  equal  to  the  mean 
difference  of  the  SNR  of  the  fused  contacts  and  the  SNR  of  the  original  contacts  that  are 
fused.  It  serves  to  give  non-fused  contacts  a  lower  detection  probability  relative  to  fused 
contacts.  (Note  that  there  is  no  physical  justification  for  this,  it  is  merely  a  heuristic  way 
to  try  and  improve  the  overall  sonar  performance  after  fusion). 

Tabic  4. 1  Different  options  to  define  true  detections  and  false  alarms  after  fusion.  TjT,  denotes  a  fused 

contact  that  originates  from  two  identical  true  targets.  T,  denotes  a  true  target  that  does  not  fuse. 
F  denotes  a  false  alarm  that  does  not  fuse.  T,F  denotes  a  fused  contact  that  originates  from  a 
true  target  and  a  false  alarm,  171^  denotes  a  fused  contact  that  originates  from  two  non-identical 
true  targets.  FF  denotes  a  fused  contact  that  originates  from  two  false  alarms. 


Option 

True  targets 

False  alarms 

Amplitude  threshold 

data  combining 

1 

Ti 

F 

T 

data  fusing 

2 

T|T,  T, 

FF  .  F  .  TtF  ,  T,Tk 

T 

3 

T|Ti  ,T|,T)F ,  TfT* 

FF.F 

T 

4 

TiTi 

FF  ,  T,F  .  TiTk 

T 

5 

TiT,  TiF ,  TiTk 

FF 

T 

12 

TiT,  T, 

FF  ,  F  ,  T,F  ,  T.Tk 

T  +  AT 

13 

T,T,  Tj.TiF ,  T|T* 

FF.F 

T  +  AT 

14 

TtTi 

FF  ,  T,F  .  T|Tk 

T  +  \T 

15 

TiT,  T,F  ,  T|Tk 

FF 

T  +  AT 

4.2.6  Pres en  ta  ti on  of  res u l ts 

For  the  single  sonars  1  and  2,  and  for  the  combined  data  sets  both  with  and  without 
fusion,  the  ROC  curves  are  generated  and  shown  by  plotting  pdi  (Ar)  versus  pjiLi  ( Ar )* 
and  pdj  (At)  versus  p/aj  (Ar). 

For  the  combined  data  set  after  fusion,  the  statistics  of  the  fusion  process  are  shown,  i.e. 
the  fraction  of  true  detections  and  false  alarms  that  are  fused.  Furthermore,  the  position 
parameters  of  the  true  detections  are  given,  both  before  and  after  fusion,  as  this  gives 
information  on  the  improvement  in  target-position  accuracy  due  to  the  fusion  process. 

4.3  Parameters  of  simulation  scenario 

For  the  simulations  the  following  parameter  values  are  chosen: 

*  Number  of  pings:  Npmg s  =  1000 

•  Sample  frequency:  100  Hz 
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Ping-repetition  rate: 
Sampled  beams: 

Sound  speed: 

St.dev.  sound  speed: 
Position  sonar  1 : 

Position  sonar  2: 

St.dev.  sonar  positions: 
St.dev.  receiver  heading: 
St.dev.  bearing  estimate: 
Time  error: 

True  position  targets: 
Mean  target  SNR: 


Or  =  2  m/s 
(xs,i ,  y,j)  =  (0 . 0)  m 
,  Vjj)  =  (1000 . 0)  m 
(o.v,,i ,  ay,,, ,  a.vVv,,- )  =  (20 . 20 , 0)  m 
<50i  =  0.5  degree 
Oq>,  =  0.5  degree 
CT/r  =0.001  S 

random  in  observed  area  around  both  sonars 
10  dB 


-90  -  90  degrees,  in  steps  of  2  degrees 
c  =  1 500  m/s 


10s 


Due  to  the  non-zero  beamspacing,  and  the  low  sample  frequency  (to  save  computation 
time),  identified  target  detections  will  in  general  not  be  exactly  at  the  true  position  of 
the  target,  but  can  be  offset  by  up  to  1 .0  degree  and  about  7.5  m. 


4.4  Single-sonar  results 


4.4. 1  Checks  on  simulations 


Before  going  into  the  possible  benefits  of  data  fusion,  we  check  the  data  generation  and 

the  method  to  measure  the  ROC  curves  by  calculating  these  curves  for 

sonars  1  and  2  separately,  and  comparing  these  to  the  theoretically  expected  ROC 

curves. 

figure  4.1  shows  the  generated  'observed'  target  positions  taken  from  bivariate  normal 
distributions,  using  the  true  target  position  and  the  position  covariances  tor  both  sonars. 
By  comparing  both  panels  in  this  Figure  it  is  clear  that  the  generated  target  positions  are 
oriented  differently  due  to  the  different  shape  and  size  of  the  covariance  ellipses  of  both 
sonars.  In  this  particular  example,  where  only  one  true  target  is  inserted,  the  positions 
are  generated  using  the  following  covariance  ellipses  around  the  true  target  position: 

•  sonar  I: 

a.v, 33.8  m 
Oy,,/:  30.1  m 

•  sonar  2: 
cx,,2:  23, 1  m 
<jy,,2-  24.3  m 

The  measured  standard  deviations  of  the  'observed'  positions  are: 

•  sonar  1: 

.E-direction:  35.0  m 
v-direction:  29.1  m 

•  sonar  2: 

.e-direction:  23.7  m 
y-direction:  24.8  m 

The  measured  standard  deviations  are  on  average  slightly  larger  than  the  standard 
deviations  of  the  generated  positions,  This  is  because  the  ‘observed’  position  of  a  target 
is  the  centre  of  the  (range,  bearing)-cell  in  which  the  generated  position  lies.  As  the 
generated  position  is  usually  not  exactly  in  the  centre  of  such  a  cell,  this  introduces  a 
(small)  extra  scatter. 
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Figure  4.2  shows  an  example  of  the  generated  target  positions  if  10  targets  are  inserted 
in  the  simulations.  This  example  serves  to  illustrate  the  fact  that  the  size  and  shape  of 
the  covariance  ellipses  depend  on  the  position  of  the  target  itself,  as  well  as  on  the 
relative  positions  of  the  source  and  receiver  sonars.  The  association  probability  will  also 
depend  on  the  position  of  the  two  contacts, 


4St  4a  «n  -«oq  -wo  -no 

X  |ffi[  Jt  jm| 

Figure  4. 1  Example  of  generated  erroneous  target  positions  for  sonar  1  (positioned  at  (x,y)  =  (0,Q)  m. 

assumed  mono  static,  left-hand  panel)  and  sonar  2  (positioned  ai  (xTy)  =  { 1000,0)  m,  assumed 
histatic.  right-hand  panel).  One  target  is  inserted  in  the  simulations.  Its  true  position  is  shown  in 
black. 


ttgure  4.2  Example  of  generated  erroneous  target  positions  for  sonar  1  (positioned  ai  (x,  y  )  =  (0,0)  mt 
assumed  monostatic,  top  panel)  and  sonar  2  (positioned  ai  (x,  y)  -  ( 1000, 0)  m,  assumed 
bistatic,  bottom  panel).  Ten  largets  are  inserted  in  the  simulations.  Their  true  positions  are 
shown  in  black. 
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De  tec  tion  pe  rfo  rman  ce 

The  detection  results  are  shown  in  Figure  4 .3  to  Figure  4,5,  Figure  4.3  shows  the 
probabilities  of  detection,  for  both  sonars,  versus  amplitude  threshold.  The  measured 
values  correspond  very  well  with  theory  (black  dashed  curve).  The  theoretical  curves 
are  given  by 


f 

pm  =  exP 


[Q'V'10  > 

2(l  +  .mr)  j 


(4.1) 


(see  e.g.  Equation  [9.4-  I4t>]  of  [7]).  Figure  4.4  shows  the  probabilities  of  false  alarm 
versus  amplitude  threshold.  Again,  the  measured  values  correspond  very  well  with 
theoretical  ones  (black  dashed  curve),  given  by 


pf„  =exp 


f 


jqAt/LO  ^ 


(4.2) 


Figure  4.5  show  s  the  ROC  curves,  i.e.  the  probability  of  detection  versus  the  probability 
of  false  alarm.  These  also  agree  very  w'cll  with  the  theoretical  curves  (black  dashed 
lines),  given  by 


_  ol/(l+w) 


P  ~  p 

1  1  fn 


(4.3) 


Note  that  the  'detection  threshold'  DT.  i.e.  the  signal-to-noise  ratio  for  P.j  =  0.5.  can  be 
derived  from  the  amplitude  threshold  A/- according  to 


(  10 


DT  =  1 0  log 

Or 


log? 
( 10  log  4 


jO'Wio  _  | 


A,  =1010log  -i^i(lOOT,10  +  l) 

1 1°  log  e 


(4.4) 


(4.5) 


giving  At  =  1 1.83  dB  if  DT  =  10. 
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ligure  4.3  Simulated  probability  of  detection  for  sonar  l  (blue)  and  2  (red),  for  a  mean  target  signabio- 
noise  ratio  of  10  dB  (relative  to  both  sonars).  The  black  dashed  lines  show  the  theoretical 
curves  for  different  assumed  SNRs  (labels  near  curves,  in  dB)  and  Rayleigh  statistics. 


Figure  4.4  Simulated  probability  of  false  alarm  for  sonar  I  (blue)  and  2  (red).  The  black  dashed  line  shows 
the  theoretical  curve  for  Rayleigh  statistics. 
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E'igure  4,5  Simulated  ROC  curves  for  sonar  1  (blue)  and  2  (red),  for  a  target  signal- to- noise  ratio  of  10  dll 
(relative  to  both  sonars),  The  black  dashed  lines  show  the  theoretical  curves  for  different 
assumed  SNRs  (labels  near  curves,  in  UB)  anti  Rayleigh  statistics. 


4.5  Statistics  of  fusions 

Before  showing  the  detection  performance  after  fusion,  we  will  first  illustrate  the  fusion 
process  itself  by  looking  at  its  statistics.  Figure  4.6  to  Figure  4.1 1  show  six  examples  of 
the  fraction  of  detections  that  fuses,  versus  threshold.  The  difference  between  the 
examples  is  the  following  parameter  values: 

•  Figure  4,6:  target  SNR:  10  dB.  multiplication  factor  for  position-error  ellipses 
(Ftuipse):  1,  number  of  targets:  !. 

•  Figure  4,7:  target  SNR:  10  dB,  I.  number  of  targets:  10. 

•  Figure  4.8:  target  SNR:  10  dB,  Ffmp!t /.  I.  number  of  targets:  30. 

■  Figure  4.9:  target  SNR:  10  dB.  multiplication  factor  for  position-error  ellipses:  2, 
number  of  targets :  10. 

•  Figure  4. 10:  target  SNR:  20  dB.  Ffuipsf:  1 .  number  of  targets:  10. 

■  Figure  4. 1 1 :  target  SNR:  20  dB.  F^iipse'  2.  number  of  targets:  10. 

In  all  Figures,  we  have  made  a  subdivision  into  five  classes: 

•  Fraction  of  false  alarms  that  fuses  with  other  false  alarms  (blue  line). 

•  Fraction  of  false  alarms  that  fuses  with  true  detections  (red  line). 

•  Fraction  of  true  detections  that  fuses  with  false  alarms  (green  line). 

•  Fraction  of  true  detections  that  fuses  with  true  detections  of  other  targets 

(e.g.  a  true  detection  of  target  2  fuses  with  a  true  detection  of  target  5.  black  line). 

•  Fraction  of  true  detections  that  fuses  with  true  detections  of  the  same  target 
(magenta  line). 

From  all  Figures  the  following  observations  can  be  made: 

•  The  fraction  of  false  alarms  that  fuses  with  other  false  alarms  (blue  lines)  decreases 
strongly  with  threshold.  This  is  because  the  number  of  false  alarms  decreases 
strongly  with  threshold,  and  thus  their  mutual  separation  increases  strongly, 
resulting  in  fewer  fusions. 
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*  The  fraction  of  false  alarms  that  fuses  with  true  detections  (red  lines)  is  small,  at 
most  1% .  It  is  smaller  than  the  fraction  of  false  alarms  that  fuses  with  other  false 
alarms.  This  is,  of  course,  as  there  are  fewer  true  detections  than  false  alarms  to  fuse 
with. 

*  The  fraction  of  true  detections  that  fuses  with  false  alarms  (green  lines)  is  about 
equal  to  the  fraction  of  false  alarms  that  fuses  with  other  false  alarms  (blue  lines). 
This  is  logical:  these  fusions  are  spurious  fusions  and  their  fraction  is  dominated  by 
the  density  of  false  alarms, 

(Note  that  this  fraction  is  in  principle  the  same  as  the  previous  one  [fraction  of  false 
alarms  that  fuses  with  true  detections],  only  scaled  with  the  total  number  of  true 
detections  instead  of  with  the  total  number  of  false  alarms.) 

•  The  fraction  of  true  detections  that  fuses  with  other  (but  different)  true  detections 
(black  lines)  is  very  small,  <1%.  This  i$  good  news:  it  show's  that,  even  if  multiple 
targets  are  present  in  the  data,  fusions  between  target  i  in  the  first  data  set  with 
target  j  *  i  in  the  second  data  set  are  expected  to  occur  not  too  often. 

•  The  fraction  of  true  detections  that  fuses  with  their  counterpart  in  the  other  data  set 
(magenta  line)  decreases  with  threshold.  The  reason  for  the  decrease  is  that  it  is  less 
likely  that  both  sonars  will  'detect*  the  target  for  higher  threshold,  and  it  is  thus  less 
likely  that  two  observations  of  the  same  target  will  fuse. 

In  Figure  4.7  (mean  target  SNR  =  10  dB,  10  targets)  the  fraction  of  true  detections  that 
fuses  with  their  counterpart  in  the  other  data  set  is  quite  low,  about  15%,  even  at  low 
threshold.  Apparently,  the  distance  between  the  two  observations  of  one  target  is  too 
large  to  allow  for  a  fusion  fraction  that  is,  say,  more  than  20%, 

One  'solution*  for  this  is  to  increase  the  size  of  the  position-error  ellipses  for  fusion 
with  respect  to  that  for  the  generation  of  the  position  errors  (see  Section  4.2.3), 

This  is  illustrated  in  Figure  4.9,  where  the  size  of  the  position -error  ellipses  for  the 
fusion  process  is  twice  as  large  as  that  for  the  generation  of  the  erroneous  target 
positions.  The  fraction  of  true  detections  that  fuses  with  their  counterpart  in  the  other 
data  set  has  increased  to  25%  at  low  thresholds. 

The  same  observation  can  be  made  by  comparing  Figure  4.10  (multiplication  factor  for 
position-error  ellipses:  1 )  and  Figure  4. 1 1  (multiplication  factor  tor  position-error 
ellipses:  2).  In  both  Figures  the  target  SNR  is  20  dB.  and  10  targets  are  inserted  in  the 
simulations.  The  maximum  fraction  of  true  fusions  increases  from  30%  (Figure  4.10) 
to  50%  (Figure  4.  II ). 

From  Figure  4.6  and  Figure  4,8,  which  both  have  a  target  SNR  of  10  dB  and  a 
multiplication  factor  for  position-error  ellipses  equal  to  1,  it  can  be  observed  that  the 
number  of  targets  in  the  simulation  (1  respectively  30)  does  not  change  significantly  the 
fraction  of  true  fusions,  although  this  fraction  increases  slightly  with  the  number  of 
targets.  It  is  to  be  expected  that,  only  if  the  number  of  targets  is  similar  to  the  number  of 
false  alarms,  the  number  of  spurious  fusions  betw  een  true  targets  and  false  alarms  w  ill 
increase,  and  the  number  of  true  fusions  will  thus  decrease. 
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Figure  4.6  Fraction  of  detections  that  fuses,  versus  threshold.  A  target  SNR  of  10  dB  is  used,  and  I  target 
is  inserted  in  the  simulations.  The  size  of  the  position -error  ellipses  for  the  fusion  process  is 
equal  to  that  for  the  generation  of  the  erroneous  target  positions.  Sir  lid  and  dashed  lines  denote 
the  mean  values  and  standard  deviations  over  all  pings. 


Amplitude  threshold  AJ  [dB] 


Figure  4.7  Fraction  of  detections  that  fuses,  versus  threshold.  A  target  SNR  of  1 0  dB  is  used,  anti  Ml 
targets  are  inserted  in  the  simulations.  The  size  of  the  position-error  ellipses  for  the  fusion 
process  is  equal  to  that  for  the  generation  of  the  erroneous  target  positions.  Solid  and  dashed 
Lines  denote  the  mean  values  and  standard  deviations  over  ail  pings. 
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Figure  4.8  Fraction  of  detections  that  fuses,  versus  threshold  A  target  SNR  of  10  dR  is  used  and  30 
targets  are  inserted  in  the  simulations.  The  si/e  of  the  posit  ton-error  ellipses  for  the  fusion 
process  is  equal  to  that  for  the  generation  of  the  erroneous  target  posit  ions.  Solid  and  dashed 
lines  denote  ihe  mean  values  and  standard  deviations  over  all  pings. 


Figure  4.9  Fraction  of  detections  that  fuses,  versus  threshold,  A  target  SNR  of  10  dR  is  used,  and  10 
targets  are  inserted  in  the  simulations.  The  size  of  the  no s i t ion -e rro r  e  1 1  i pse s  for  the  fusion 
process  is  two  times  lamer  than  for  the  generation  of  the  erroneous  target  positions.  Solid  and 
dashed  lines  denote  the  mean  values  and  standard  deviations  over  all  pings. 
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figure  4  JO  Fraction  of  detections  ihut  luses.  versus  threshold.  A  target  SNR  of  20  dB  is  used,  and  10 
targets  are  inserted  in  the  simulations.  The  size  of  the  posit  ion -error  ellipses  tor  the  fusion 
process  is  equal  to  that  for  the  generation  of  the  erroneous  target  positions.  Solid  and  dashed 
lines  denote  the  mean  values  and  standard  deviations  over  all  pings. 
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Figure  4, 1 1  Fraction  of  detections  that  fuses,  versus  threshold.  A  target  SNR  of  20  dB  is  used,  and  1 0 
targets  are  inserted  in  the  simulations.  The  si/.e  of  the  position-error  ellipses  for  the  fusion 
process  is  two  times  larger  than  that  for  the  generation  of  the  erroneous  target  positions.  Solid 
and  dashed  lines  denote  the  mean  values  and  standard  deviations  over  all  pings. 


4.6  Increase  in  amplitude  threshold  after  fusion 

As  was  mentioned  in  Section  4,2.5,  we  allow-  the  option  to  increase  the  amplitude 
threshold  based  on  the  statistics  of  the  fusion  process.  That  is.  when  two  contacts  are 
fused  the  SNR  of  the  fused  contact  is  set  equal  to  the  maximum  of  the  two  original 
contacts.  This  means  that  the  mean  SNR  assigned  to  all  fused  contacts  will  he  higher 
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than  the  mean  original  SNR  of  all  contacts  that  are  fused.  This  increase  is  denoted  by 
AAy.  In  order  to  give  higher  confidence  to  the  fused  contacts  in  the  detection  process, 
we  can  increase  the  amplitude  threshold  by  AAy,  so  that  fewer  non-fused  contacts  will 
pass  the  new  threshold  Ar  +AAy, 

Figure  4. 1 2  shows  AAy  versus  A r  for  the  case  where  there  are  10  targets,  each  with  a 
mean  SNR  of  10  dB.  It  is  clear  from  this  Figure  that  at  low  thresholds  ( Ar  <  15  dB) 
the  increase  in  threshold  is  quite  large,  sometimes  more  than  1  dB.  This  means  that 
instead  of  applying  the  nominal  threshold  of  e.g.  11  dB,  a  threshold  of  1 1 .5  dB  is 
applied.  This  will  only  occasionally  affect  the  number  of  fused  contacts  that  are 
detected,  but  it  will  decrease  the  number  of  non-fused  contacts  for  that  threshold. 

It  should  be  remarked  that; 

•  A Ar  is  significantly  larger  for  the  fusion  of  true  contacts  than  for  the  fusion  of  false 
contacts.  This  is  because  the  mean  SNR  of  the  true  contacts  is  10  dB  larger  than  the 
mean  SNR  of  false  alarms.  With  this  higher  mean  SNR.  the  spread  will  also  be 
larger. 

•  A A/-  decreases  with  Ay  for  both  categories.  This  is  because  the  difference  in  SNR 
between  both  original  contacts  will  decrease  with  Ar-  Both  contacts  have  to  have  an 
SNR  larger  than  Ay,  and  the  higher  Ay,  the  smaller  the  difference  in  both  SNRs. 

The  mean  AAy  (blue  line)  increases  between  Ay  =  10  and  Ay  =  13  dB.  because  the 
relative  ratio  of  true  -  true  and  false  -  false  fusions  increases  with  Ay. 


Figure  4.12  Increase  in  amplitude  threshold,  calculated  from  ihe  mean  difference  in  SNR  he  lure  and  alter 
fusion  (see  text  above  Figure  for  details).  Solid  lines:  moan  increase  in  threshold  over  all  pings, 
dashed  lines:  standard  deviation.  Red  lines:  only  using  false-alarm  -  false-alarm  fusions,  green 
line:  only  using  true -detection  -  true -detection  fusions,  blue  line:  using  all  fusions. 


47  False-aiarm  probability  after  fusion 

Figure  4, 13  shows  the  false-alarm  probabilities  after  fusing  the  two  data  sets. 

The  number  of  targets  inserted  in  the  simulations  is  10.  and  their  mean  SNR  is 
10  dB.  The  size  of  the  position-error  ellipses  for  the  fusion  process  is  equal  to  that  for 
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the  generation  of  the  erroneous  target  positions.  The  different  options  correspond  to 
those  in  Table  4.1. 

The  following  observations  can  be  made  from  Figure  4.13: 

•  The  faise-alarm  probabilities  of  the  combined  data  set  and  for  options 

2.12,  and  3.  tire  a  factor  two  larger  than  that  of  the  two  separate  sonars  (as  shown  in 
Figure  4.4).  This  is  because  the  total  number  of  false  alarms  is  not  scaled  with  the 
total  number  of  faise-alarm  opportunities  in  both  data  sets,  but  with  the  total 
number  of  faise-alarm  opportunities  in  one  data  set. 

•  The  fall-off  of  the  false  alarm  probability  is  the  same  as  the  theoretical  curve  (black 
dashed  line)  that  is  based  on  Rayleigh  statistics.  There  is  no  reason  for  the  fused 
contacts  to  follow  this  curve,  but  as  there  are  more  non-fused  than  fused  contacts, 
the  non-fused  contacts  dominate  the  statistics. 

•  The  faise-alarm  probabilities  for  options  4.  14.  5  and  15  (i.e.,  only  considering 
fused  contacts)  are  a  factor  of  10  -  30  smaller  than  those  for  options 

2.  12.  3  and  13  (which  consider  all  contacts  after  fusion). 

•  The  faise-alarm  probabilities  for  options  12  and  13  (i.e..  above  A/  +A Ar )  are 
roughly  a  factor  of  two  smaller  than  those  for  option  2  and  3  (i.e.,  above  Aft, 

The  number  of  false  alarms  is  thus  decreased  significantly  by  raising  the  amplitude 
threshold  by  A A/- 

•  The  faise-alarm  probabilities  for  options  14  and  15  (i.e.,  above  At+&At  )  are  only 
slightly  smaller  than  those  for  option  4  and  5  (i.e.,  above  4/-).  That  is.  the  number  of 
false  alarms  is  not  decreased  significantly  by  raising  the  amplitude  threshold  if  only 
fused  contacts  are  considered. 

For  the  influence  of  the  number  of  targets  inserted  in  the  simulations,  their  mean  SNR. 
and  the  size  of  the  position-error  ellipses  on  the  faise-alarm  probabilities, 
see  Section  4.10. 


Figure  4.13  False-ilarm  probability  for  sonars  1  and  2  combined  (blue  solid  line)  and  sonars  1  and  2  fused 
(other  lines).  The  different  options  are  explained  in  Table  4. 1  (solid  lines,  kURf  fusion,  dashed 
lines:  '  AND'  fusion).  The  results  for  options  2  and  3  overlap,  as  well  as  those  for  options  12 
and  13.  For  options  4  and  5  and  for  options  14  and  15.  The  black  dashed  line  shows  the 
theoretical  curve  for  Rayleigh  statistics. 
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4,8  Probability  of  detection  after  fusion 

Figure  4.14  shows  the  probability  of  detection  after  fusing  the  two  data  sets.  In  both 
examples,  the  number  of  targets  inserted  in  the  simulations  is  10,  and  their  mean  SNR  is 
10  dB,  The  size  of  the  position-error  ellipses  for  the  fusion  process  is  equal  to  that  for 
the  generation  of  the  erroneous  target  positions.  The  different  options  correspond  to 
those  in  Table  4.1,  The  following  observations  can  be  made  from  Figure  4, 14; 

*  The  probability  of  detection  of  the  combined  data  set  (blue  solid  line)  is  nearly 
equal  to  the  theoretical  value  after  ‘OR*  fusion, 

(4.6) 

where  Pi  is  the  probability  of  detection  of  sonar  i  separately  (black  dotted  lines). 
This  is  expected  from  theory,  for  two  independent  observations. 

*  The  probabilities  of  detection  for  options  2,  and  3  are  almost  equal  to  that  of  the 
combined  data  set.  Those  for  options  12  and  13  are  slightly  smaller,  about  0,02, 

*  The  probabilities  of  detection  for  options  4,  14,  5  and  15  (i.e..  only  taking  into 
account  fused  contacts)  are  much  smaller  than  those  of  the  combined  data  set,  and  a 
factor  2-3  smaller  than  the  theoretical  curves  (black  dashed  lines).  The  latter  is 
due  to  the  small  number  of  fusions  that  occurs  (see  Section  4.5).  which  results  from 
the  position  errors  on  the  single-sonar  contacts.  Without  such  position  errors,  more 
fusions  would  occur  and  the  results  would  follow  the  theoretical  lines  much  better. 

*  Including  the  increased  threshold  At+&At  instead  of  At  decreases  the  probability  of 
detection  by  about  0.02. 

For  the  influence  of  the  number  of  targets  inserted  in  the  simulations,  their  mean  SNR. 
and  the  size  of  the  position-error  ellipses  on  the  probability  of  detection, 
see  Section  4.10. 


Figure  4,14  Probability  of  detection  for  sonars  l  and  2  combined  (blue  solid  line)  and  sonars  1  and  2  fused. 

The  different  opiums  are  explained  in  Table  4. 1  (solid  lines:  lORb  fusion,  dashed  lines:  'AND1 
fusion).  Black  dotted  lines:  theoretical  curves  for  'OR1  fusion  and  different  assumed  SNRs, 
Black  dashed  lines:  theoretical  curves  for  ‘AND’  fusion  and  different  assumed  SNRs. 
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4.9  ROC  curves  after  fusion 

Figure  4.15  shows  the  ROC  curves  after  fusing  the  two  data  sets.  In  both  examples,  the 
number  of  targets  inserted  in  the  simulations  is  10.  and  their  mean  SNR  is 
10  dB.  The  size  of  the  position-error  ellipses  for  the  fusion  process  is  equal  to  that  for 
the  generation  of  the  erroneous  target  positions.  The  different  options  correspond  to 
those  in  Table  4.1. 

The  following  observations  can  be  made  from  Figure  4.15: 

•  The  operational  SNR'  of  the  combined  data  set  (blue  solid  line)  is  significantly 
higher  than  that  of  the  separate  data  sets  (see  Figure  4.16).  This  is  a  considerable 
performance  improvement. 

■  The  operational  SNRs  of  the  fused  data  sets  with  options  2.  12.  3  and  13  is  also 
t  -  5  dB  higher  than  that  of  the  separate  data  sets,  and  similar  to  that  of  the 
combined  data  set. 

•  The  ROC  curves  for  the  combined  data  set  and  of  the  fused  data  sets  with 
options  2.  12.  3  and  13  closely  follow  the  theoretical  ROC  curves  for  the  'OR' 
fusion  (black  dotted  lines).  The  theoretical  ROC  curves  are  obtained  by  plotting 
l  -( I  -PdnY  versus  1  -( 1  -P)aY.  where  P,ui  and  P/a  are  given  by  Equations  (4. 1 )  and 
(4.2),  respectively. 

•  The  operational  SNRs  of  the  fused  data  sets  with  options  4,  14,  5  and  15  are  much 
below  that  of  the  combined  data  set.  The  main  reason  is  the  small  number  of  fusions 
that  occurs  (see  Section  4.5),  which  is  due  to  the  position  errors  on  the  single-sonar 
contacts.  Without  such  position  errors,  more  fusions  would  occur. 

•  The  ROC  curves  of  the  fused  data  sets  with  options  4,  14,  5  and  1 5  are  reasonably 
close  to  the  theoretical  ROC  curves  for  the  'AND'  fusion  (black  dashed  lines). 

The  theoretical  ROC  curves  are  obtained  by  plotting  P,^2  versus  P/„2.  where  P,{„ 
and  P/„  are  given  by  Equations  (4.1)  and  (4.2),  respectively. 

For  the  influence  of  the  number  of  targets  inserted  in  the  simulations,  their  mean  SNR. 
and  the  size  of  the  position-error  ellipses  on  the  ROC  curves,  see  Section  4. 10. 


J  The  operational  SNR  is  usually  defined  as  the  SNR  of  the  theoretical  ROC  curves  that  corresponds  best 
to  the  data 
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I  igure  4.15  ROC  curves  tor  sonars  I  and  2  combined  (blue  solid  line)  and  sonars  1  and  2  fused,  The 
different  options  are  explained  in  Table  4.1  (solid  lines:  “OR*  fusion,  dashed  lines;  “AND1 
fusion).  10  targets  with  mean  SNR  =  10  dll  are  inserted.  Black  dotted  lines;  theoretical  curves 
for  “OR'  fusion  and  different  assumed  SNRs.  Black  dashed  lines:  theoretical  curves  for  'AND' 
fusion  and  different  assumed  SNRs. 


figure  4. 16  Comparison  between  single-sonar  ROC  curves  (black  lines)  and  ROC  curves  after  combining 
or  fusion  the  contacts  tall  other  colors).  10  targets  with  mean  SNR  -  10  dB  are  inserted. 


Summary 

The  above  results  can  be  summarised  as  follows: 

Combining  or  ‘OR1  fusing  two  data  sets,  with  a  Rayleigh  distribution  tor  the  target 
amplitudes,  gives  a  considerable  improvement  of  performance  over  using  a  single 
sonar.  In  hindsight  this  is  to  be  expected  and  it  can  be  derived  from  a  proper  theoretical 
analysis. 

Using  only  the  fused  contacts  ('AND’  fusion),  the  false-alami  probability  is  very  small 
and  reduces  to  almost  zero  for  amplitude  thresholds  >  12  dB.  However,  the  probability 
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of  detection  also  reduces  from  about  0.6  to  0.  L  This  may  be  due  to  the  low  mean  SNR 
of  the  target  ( 10  dB  in  this  case). 

Using  both  the  fused  and  non-fused  contacts  after  fusion  ('OR'  fusion),  the  probability 
of  detection  and  false  alarm  are  roughly  equal  to  those  if  the  data  sets  are  Just' 
combined  (i.e.,  with  no  fusion  applied).  ’OR"  fusion  thus  does  not  seem  to  give  much 
benefit  over  combining  the  two  data  sets.  However,  in  Section  4,1 1  we  w  ilt  show-  that  in 
some  cases  it  may  still  be  desirable  to  apply  the  fusion. 

Before  showing  the  benefits  of  fusion  in  Section  4.1  U  the  next  Section  show  s  the  ROC 
curves  for  different  values  of  the  input  parameters  of  the  simulation, 

4.10  ROC  curves  for  different  simulation  parameters 

To  study  the  dependence  of  the  fusion  on  the  input  parameters  of  the  simulation,  we 
have  varied  the  following  parameters: 

*  Number  of  targets:  1,  3,  10,  30. 

*  Mean  target  SNR:  10,  15,  20  dB, 

*  Size  of  error  ellipse  used  tor  fusion  (see  Section  4,2.3);  1,  2,  3  times  larger  than  size 
of  error  ellipse  used  to  generate  the  erroneous  target  positions. 

Each  of  these  three  variations  is  treated  in  one  of  the  following  three  Sections. 

4.1QJ  Number  of  targets 

Figure  4. 1 7  shows  the  ROC  curves  for  different  numbers  of  targets  ( 1  and  30). 

The  mean  target  SNR  is  set  to  10  dB,  and  the  size  of  the  error  ellipses  is  equal  to  that 
used  for  generating  the  erroneous  target  positions.  It  is  found  that  the  number  of  targets 
inserted  in  the  simulations  does  not  significantly  influence  the  false-alarm  probabilities 
or  the  probability  of  detection,  as  is  expected.  A  larger  number  of  targets  only  leads  to 
better  statistics,  but  does  not  significantly  increase  the  number  of  fusions  between 
contacts  that  originate  from  different  targets. 


Figure  4. 17  ROC  curves  for  sonars  l  and  2  combined  iblue  solid  line)  and  sonars  I  and  2  fused.  The 

different  options  are  explained  in  Table  4.  L  Left:  I  target.  Right:  30  targets.  Black  dotted  lines: 
theoretical  curves  for  ‘OR*  fusion  and  different  assumed  SNRs.  Black  dashed  lines:  theoretical 
curves  for  ‘AND'  fusion  and  different  assumed  SNRs. 


4.10.2  Target  SNR 

It  is  found  that  the  target  SNR  does  influence  the  false-alarm  probabilities  if  the  mean 
SNR  is  quite  high.  e.g.  20  dB  (see  Figure  4, 18).  The  reason  for  this  is  the  fusion  of 
target  contacts  with  false  alarms  (TjF  option  in  Table  4. 1 ),  In  fact,  only  about  one  such 
fusion  in  100  simulations  is  enough  to  account  for  this. 
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The  target  SNR  does  influence  the  probability  of  detection  quite  strongly,  as  expected 
(see  Figure  4.19).  In  this  Figure  the  number  of  targets  is  10.  and  the  size  of  the  error 
ellipses  is  equal  to  that  used  for  generating  the  erroneous  target  positions.  The  curves 
for  options  2,  12.  3,  and  13  follow  nicely  the  theoretical  curves  for  ‘OR'  fusion. 

The  curves  for  options  4,  14.  5  and  15  are  much  lower  than  the  theoretical  curves  for 
'AND'  fusion,  due  to  the  position  errors  in  the  original  contacts. 

Figure  4.20  shows  the  resulting  ROC  curves  for  different  mean  target  SNRs 
(10  and  20  dB).  The  ROC  curves  for  all  options  improve  if  the  target  SNR  increases. 
For  higher  target  SNR.  the  ROC  curves  deviate  more  and  more  from  the  theoretical 
curves  and  some  of  the  curves  show  peculiar  patterns,  The  reason  for  this  is  the 
relatively  large  false-alarm  probabilities  at  high  thresholds  (see  Figure  4. 18). 


Figure  4. 1 8  False-alarm  probabilities  for  sonars  1  and  2  combined  (blue  solid  line)  and  sonars  I  and  2 
fused  The  different  options  are  explained  in  Table  4.1,  Left:  mean  target  SNR  =  10  dll. 
Right;  mean  target  SNR  =  20  dB.  Black  dotted  line:  theoretical  curve  for  Rayleigh  statistics. 


Figure  4.  ]  0  Probability  of  detection  for  sonars  1  and  2  combined  (blue  solid  line)  and  sonars  l  and  2  fused. 

The  different  options  are  explained  in  Table  4.1. 1  .eft :  mean  target  SNR  =  10  dB.  Right:  mean 
target  SNR  =  20  dB,  Black  dotted  lines:  theoretical  curves  for  ’OR*  fusion  and  different 
assumed  SNRs  (see  labels).  Black  dashed  lines:  theoretical  curves  for  ‘AND ‘  fusion  and 
different  assumed  SNRs  (see  labels). 
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Figure  4.20  ROC  curves  for  sonars  l  and  2  combined  (blue  solid  line)  and  sonars  1  and  2  fused. 

The  different  options  are  explained  in  Table  4.1.  Left:  mean  target  SNR  =  10  dB,  Right:  mean 
target  SNR  =  20  dB,  Blade  dotted  lines:  theoretical  curves  for  ’OR'  fusion  and  different 
assumed  SNRs  (left:  K  -  12  dB,  right:  15  -  20  dB),  Black  dashed  lines:  theoretical  curves  for 
'AND'  fusion  and  different  assumed  SNRs  (left:  5  10  dB.  right:  10  15  dB). 


4.  10.3  Size  of  error  ellipse  used  for  fusion 

It  is  found  that  the  size  of  the  position-error  ellipses  for  the  fusion  process  does  not 
significantly  influence  the  false-alarm  probabilities.  Increasing  the  size  of  the  position- 
error  ellipses  by  a  factor  of  3  roughly  doubles  the  probability  of  detection  for 
options  4.  14.  5  and  15  (see  Figure  4.21).  In  this  Figure  the  number  of  targets  is  10,  and 
the  target  mean  SNR  is  set  to  10  dB.  Figure  4.22  shows  the  corresponding  ROC  curves 
for  different  sizes  of  the  error  ellipses  (factor  1  and  3  larger  than  error  ellipse  used  to 
generate  target  positions). 

The  ROC  curves  only  improve  for  options  4,  14*  5  and  15  (which  all  only  include  fused 
contacts),  and  mainly  at  lower  amplitude  thresholds.  The  reason  tor  this  is  that  it 
becomes  easier  for  contacts  to  fuse. 


Figure  4,21  Probability  of  detection  for  sonars  1  and  2  combined  (blue  solid  line)  and  sonars  I  and  2  fused. 

The  different  options  are  explained  in  Table  4. 1 .  The  size  of  the  error  ellipse  used  for  fusion  is 
a  factor  l  (left)  and  3  (right)  larger  than  the  size  of  ihe  error  ellipse  used  to  generate  the 
erroneous  target  positions.  Black  dotted  lines:  theoretical  curves  lor  ‘OR*  fusion  and  different 
assumed  SNRs,  Black  dashed  lines:  theoretical  curves  for  ‘AND*  fusion  and  different  assumed 
SNRs. 
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Figure  4.22  ROC  curves  for  sonars  l  and  2  combined  (  blue  solid  line  )  and  sonars  l  and  2  fused.  The 
different  options  are  explained  in  Table  4. 1 ,  The  si/e  of  the  error  ellipse  used  for  fusion  is 
a  factor  \  (left)  and  3  (right  )  larger  than  the  size  of  the  error  ellipse  used  to  generate  the 
erroneous  target  positions.  Black  dotted  lines:  theoretical  curves  for  ’OR'  fusion  and  different 
assumed  SNRs,  Black  dashed  lines:  theoretical  curves  for  *  AND'  fusion  and  different  assumed 
SNRs. 


4. 10.4  Target  amplitude  distribution 

So  far,  a  Rayleigh  distribution  was  assumed  for  the  background  and  for  the  target 
amplitude.  This  basically  means  that  the  fluctuations  in  the  amplitude  of  both  are 
random  in  nature.  An  alternative  distribution,  used  for  radar  targets  containing  a  strong 
highlight,  is  known  as  the  one-dominam-plus- Rayleigh  (ID+R)  distribution  [7], 

We  have  evaluated  the  benefits  of  fusion  with  this  ID+R  distribution  as  well,  and 
compared  it  to  the  benefits  of  a  Rayleigh  distribution,  see  Figure  4,23,  Although  the 
ROC  curves  for  both  distributions  are  different,  and  the  quantitative  change  due  to 
fusion  is  different  for  both  statistics,  the  overall  improvement  or  decrease  between  the 
single-sonar  performance  (black  lines)  and  that  after  fusion  (lines  in  other  colors)  is 
similar  for  both  situations. 


00 
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Figure  4.23  Comparison  between  single-sonar  ROC  curves  (black  lines!  and  ROC  curves  after  combining 
or  fusion  the  contacts  (all  other  colors)*  using  10  targets  with  a  mean  SNR  of  10  dR. 

Left;  Rayleigh  distribution  of  target  amplitude. 

Right:  One-do  rmnant+  Ray  lei  gh  d  i  stri  butio  n  of  t  argei  amp  I  i  tude . 


4.1 1  Improvement  in  localisation  after  fusion 

The  target  detections  with  one  sonar  are  usually  observed  at  inaccurate  positions 
(see  also  Section  4.4),  Due  to  (small)  errors  in  the  sonar  position,  receiver  heading, 
sound  speed  and  time,  the  detections  de  viate  from  the  true  position  of  the  target 
(see  Figure  4. 1  and  Figure  4.2  for  an  illustration).  Apart  from  the  increase  in  detection 
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performance  after  fusion  (see  previous  Sections),  another  benefit  of  fusing  two 

(or  more)  data  sets  is  the  increased  position  information  that  one  obtains  for  the  targets. 

For  the  simulations  with  30  targets  inserted,  this  is  illustrated  in  Table  4.2. 

For  sonars  1  and  2 ,  'Theory'  gives  the  parameters  of  the  covariance  ellipses  according 
to  which  the  target  positions  are  generated.  'Generated  positions'  refers  to  the  positions 
where  targets  are  placed  in  the  simulations,  given  the  errors  in  some  fundamental 
parameters.  These  agree  well  with  the  theoretical  values. 

'Observed  positions'  refers  to  the  position  where  the  targets  are  observed,  given  the 
beam  spacing  and  time  sampling  of  the  sonars.  As  explained  in  Section  4.41,  these  are 
slightly  larger  than  the  theoretical  errors.  This  is  because  the  observed  position  of  a 
target  is  die  centre  of  the  (range  ,  bearmg)-eell  in  which  the  generated  position  lies. 

As  the  generated  position  is  usually  not  exactly  in  the  centre  of  such  a  cell,  this 
introduces  a  (small)  extra  scatter. 

The  following  remarks  can  be  made  from  this  Table: 

•  The  mean  offsets  of  the  generated  and  observed  positions  of  sonars  1  and  2  are 
always  small,  1  -  3  metres, 

•  Using  a  fixed  threshold  of  10  dB,  and  evaluating  fusion  options  2,  12,  3  and  13,  it 
can  be  witnessed  that  the  standard  deviations  of  the  observed  positions  after  fusion 
are  similar  to  those  before  fusion.  Most  contacts  included  in  these  cases  are 
non-fused  contacts,  which  have  the  same  position  accuracy  as  those  before  fusion, 

•  Using  a  fixed  threshold  of  10  dB,  and  evaluating  fusion  options  4  and  14  (and  to  a 
lesser  extent  options  5  and  15),  it  can  be  witnessed  that  the  standard  deviations  of 
the  observed  positions  after  fusion  are  smaller  than  those  before  fusion:  the  position 
estimate  of  contacts  that  originate  from  the  fusion  of  two  true  contacts  have  a 
position  accuracy  that  is  roughly  a  factor  1,5  times  better  than  that  of  the  original 
contacts.  This  is  an  important  observation  that  will.  e,g.  provide  better  input 
parameters  for  a  tracking  algorithm. 

•  For  options  2  and  3  (“OR*  fusion),  the  standard  deviations  of  the  observed  positions 
decrease  with  threshold  Ar ♦  especially  at  higher  thresholds,  Af>  16  dB. 

The  standard  deviations  have  roughly  reduced  to  one  third  of  their  original  value 
between  At  =  10  and  Aj ■=  20  dB.  This  is  because  at  higher  threshold  relatively  more 
fused  contacts  are  present,  that  have  a  higher  position  accuracy; 

•  For  options  4  and  5  ('AND'  fusion),  the  standard  deviations  of  the  observed 
positions  also  decreases  with  threshold  Ar*  the  decrease  being  more  profound  at 
lower  values  of  Aj  than  for  options  2  and  3. 

These  results  are  found  to  be  independent  of  the  size  of  the  error  ellipses  used  for  the 
fusion. 
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Table  4.2  Summary  of  position  emirs  of  generated  and  observed  target  positions  relative  to  true  target 
positions.  All  values  are  in  metres  and  are  mean  values  over  10  targets  and  UK)  pings.  N  is  the 
total  number  of  contacts  used  in  calculating  the  position  errors. 


x-direction 

y-direction 

Mean 

St.dev. 

Mean 

St.dev. 

Sonar  1 

Theory 

0 

40 

0 

48 

Generated  positions 

2 

41 

-2 

49 

Observed  positions 

2 

45 

1 

54 

Sonar  2 

Theory 

0 

33 

0 

39 

Generated  positions 

-1 

31 

1 

37 

Observed  positions 

-2 

35 

3 

44 

After  ‘OFT  fusion  (Ar  =  10  dB) 

Option  2  (A/  =  1116) 

0 

43 

2 

55 

Option  12  (N  =  1078} 

0 

42 

3 

56 

Option  3  (N=  1154) 

0 

42 

2 

56 

Option  13  (N=  1116} 

0 

42 

3 

56 

After  'AND'  fusion  10  dB) 

Option  4  (N=  128) 

0 

23 

1 

34 

Option  14  (W=  128) 

0 

23 

1 

34 

Option  5  (N  =  193) 

0 

30 

’4 

4B 

Option  15  (A/  -  193) 

0 

30 

-4 

48 

After  fusion  (Option  2,  'OR'} 

10  dB  (N=  1116) 

0 

42 

2 

55 

Ar=12dB(N  =  869) 

0 

40 

3 

52 

Ar=  14  dB  (N=  589) 

-1 

40 

4 

51 

Ar*  16  dB  (N  =  259) 

-1 

36 

5 

47 

At=  18  dB  (N=  85) 

0 

28 

4 

30 

At  =20  6B  (N=  24} 

-2 

15 

2 

19 

At=  22  dB  (A/  =  1) 

-1 

0 

1 

0 

After  fusion  (Option  3,  ‘OR’) 

At=  10  dB  {N=  1154} 

0 

42 

2 

56 

At=  12  dB  (N=  869) 

0 

39 

3 

51 

At=  14  dB  (N=  589) 

-1 

39 

4 

50 

At=  16  dB  {N=  259) 

-1 

35 

5 

45 

At=  18  d8  (N=  85) 

0 

25 

3 

27 

At  -  20  dB  {N  =  24) 

-1 

13 

2 

17 

Ar=  22  dB  (N  =  1) 

-1 

0 

0 

0 

After  fusion  (Option  4,  ’AND’) 

Ar=  10  dB  (N=  128) 

0 

23 

1 

34 

At=  12  dB  (N=  106) 

0 

21 

0 

32 

At  =  14  dB  (N  =46) 

2 

23 

-2 

27 

At=  16  dB  (N=  9) 

-3 

15 

4 

26 

At=  17  dB  {N=  5) 

-5 

12 

5 

21 

After  fusion  (Option  5h  ‘AND') 

At=  10  dB  (N  =  193) 

0 

30 

-4 

46 

At=  12  dB  (N=  106) 

0 

21 

0 

32 

At=  14  dB  (W=  46) 

2 

23 

-2 

27 

At=  16  dB  (N=  9) 

-3 

15 

4 

26 

At  -  17  dB  (W  =  5) 

*5 

12 

5 

21 
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